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ABSTRACT 

This study explores the effects of different assumptions and systematics on the determination of the 
local, spatially resolved star formation law. Using four star formation rate (SFR) tracers (Ha with 
azimuthally averaged extinction correction, mid-infrared 24 /im, combined Ha and mid-infrared 24 
/im, and combined far- ultraviolet and mid- infrared 24 fj.m), several fitting procedures, and different 
sampling strategies we probe the relation between SFR and molecular gas at various spatial resolutions 
(500 pc and larger) and surface densities (Sh 2 ~ 10 — 245 M© pc -2 ) within the central ~ 6.5 kpc in 
the disk of NGC 4254. We explore the effect of diffuse emission using an unsharp masking technique 
with varying kernel size. The fraction of diffuse emission, /de, thus determined is a strong inverse 
function of the size of the filtering kernel. We find that in the high surface brightness regions of 
NGC 4254 the form of the molecular gas star formation law is robustly determined and approximately 
linear (~ 0.8 — 1.1) and independent of the assumed fraction of diffuse emission and the SFR tracer 
employed. When the low surface brightness regions are included, the slope of the star formation 
law depends primarily on the assumed fraction of diffuse emission. In such case, results range from 
linear when the fraction of diffuse emission in the SFR tracer is /de ^ 30% (or when diffuse emission 
is removed in both the star formation and the molecular gas tracer), to super- linear (~ 1.4) when 
/de <; 50%. We find that the tightness of the correlation between gas and star formation varies with 
the choice of star formation tracer. The 24 /zm SFR tracer by itself shows the tightest correlation 
with the molecular gas surface density, whereas the Ha corrected for extinction using an azimuthally- 
averaged correction shows the highest dispersion. We find that for R < O.5R25 the local star formation 
efficiency is constant and similar to that observed in other large spirals, with a molecular gas depletion 
time r dep ~ 2 Gyr. 

Subject headings: galaxies: general — galaxies: individual(NGC 4254) — galaxies: spiral — galaxies: 
ISM — ISM:molecules — stars: formation 



1. INTRODUCTION 

The formation and evolution of galaxies is driven by 
the complex processes of star formation (SF) that oc- 
cur inside them. Some galaxies produce stars at very 
low rates < 0.1 M Q yr _1 , some do at modest rates ~ 1 
M© yr" 1 , while others host ongoing starbursts with SFR, 
i~ 10 — 1000 M© yr _1 . The processes responsible for con- 
verting gas into stars in various galactic environments are 
still poorly understood. Observations find that the SFR 
and the gas content in galaxies are related by, 

Ssfr = A Sg as , (1) 
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where Ssfr and E gas are the star formation rate surface 
density and the gas (atomic and molecular) surface den- 
sity, respectively; and A is the normalization constant 
representing the efficiency of the processes (Schmidt 
1959, 1963; Sanduleak 1969; Hartwick 1971; Kennicutt 
1989). For disk averaged surface densities, both normal 
star-forming and starburst galaxies follow this relation- 
ship with a power law index of N ~ 1.3 — 1.5 for total gas 
(Kennicutt 1989, 1998a, b). This relationship between 
the gas and SFR surface densities is commonly referred 
to as the Schmidt-Kennicutt SF law. This is in principle 
consistent with large scale gravitational instability being 
the major driver (Quirk 1972; Madore 1977). 

Although spatially unresolved studies of HI, CO, and 
SFR are useful for characterizing global disk properties, 
understanding the mechanisms behind the SF law re- 
quires resolved measurements. Only recently has it be- 
come possible to probe the form of the gas-SF relation- 
ship on kpc and sub-kpc scales, through the availability 
of high-resolution interferometric HI and single-dish CO 
observations and of a suite of multi- wavelength SFR trac- 
ers. Studies of the local SF law on nearby galaxies pro- 
vide substantial evidence that the molecular gas is well- 
correlated with the SFR tracers, whereas the atomic gas 
shows little or no correlation with SF activity (e.g., Wong 
& Blitz 2002; Bigiel et al.. 2008). This is a natural con- 
sequence of stars forming out of giant molecular clouds 
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(GMCs), as we observe in the local universe. Moreover, it 
has long been known that the spatial distribution of CO 
emission follows closely that of the stellar light and Ha 
(Young & Scoville 1982; Scoville k Young 1983; Solomon 
et al. 1983; Lord & Young 1990; Tacconi & Young 1990; 
Boselli et al. 1995). 

Spatially resolved SF law studies frequently reach dis- 
similar conclusions on the value of the exponent in Equa- 
tion [T] when relating molecular gas to SFR. Hereafter we 
will express the exponent as N mo i to represent the molec- 
ular gas SF law. Wong & Blitz (2002) used azimuthally 
averaged radial profiles for gas and SFR in a sample of 
seven molecule-rich spiral galaxies, finding that the best 
fit power law index for the molecular gas and SFR density 
radial profiles is N mo i ~ 0.8 — 1.4, very dependent on the 
extinction correction applied to their SFR tracer (Ha). 
Boissier et al. (2003) used CO observations along the ma- 
jor axes of sixteen disk galaxies with spatial resolution of 
— 1 — 4 kpc to carry out a similar study. They found a 
somewhat steeper exponent N mo i — 1 — 2.3 with respect 
to the molecular gas, but with a different choice for the 
extinction correction and different fitting methodologies. 
Sampling various spatial scales and surface densities of a 
sample of twenty three disk galaxies Komugi et al. (2005) 
found N roo i — 1.3 using extinction corrected Ha as the 
SFR tracer. In the most recent comprehensive study, 
Bigiel et al. (2008) analyzed a sample of 18 normal disk 
and irregular galaxies using a combination of GALEX 
far-ultraviolet (FUV, 1350-1750 A) emission and Spitzer 
24 /im to trace SFR, and CO 2 — 1 for the molecular gas 
at a spatial resolution — 750 pc. They found a best fit 
slope N mo i «1.0 for the power-law relation between Sh 2 
and Esfr- 

This spread in the value of the power law index is ob- 
served by in-depth case studies of just one galaxy. An ex- 
ample is M 51, where Kennicutt et al. (2007) studied the 
relation between gas and SFR on ~ 0.3 — 1.8 kpc scales 
sampling the distribution of emission in circular aper- 
tures centered on H II regions to find N mo i ~ 1.37— 1.56, 
depending on the spatial scale considered. Another ex- 
ample is NGC 7331 studied by Thilker et al. (2007). The 
authors find N mo i ~ 1.64 using bolometric (combining 
ultraviolet and infrared luminosity) SFR tracer at 400 
pc resolution. By comparison, Schuster et al. (2007) 
used the A = 20 cm radio-continuum as SFR tracer and 
CO J = 2 — 1 data to find a much shallower molecular 
gas power-law N mo i < 1, changing with galactocentric 
distance. Similarly, Blanc et al. (2009) studied the cen- 
tral —4 kpc of M 51 using optical spectroscopic Ha data 
at 170 pc resolution, finding a slightly sub-linear rela- 
tionship (N mo i — 0.82 ± 0.05). A very recent example 
is M 33, where Verley et al. (2010) employed a range of 
methods of data sampling, fitting techniques, and SFR 
tracers to find that the functional form of molecular gas- 
SFR relation varies (N mo i ~ 1.2 — 1.8) depending on the 
choice of SFR tracers, data sampling and fitting tech- 
niques. For the same galaxy but with a different SFR 
tracer Heyer et al. (2004) found N mol — 1.4. 

The spread in the value of the power-law index within 
and among galaxies may be intrinsic and contain valuable 
astrophysical information, or be entirely attributable to 
the different choices of gas and SFR tracers, method- 
ologies for internal extinction correction, differences in 



the CO-to-H2 conversion, or the range of spatial scales 
probed. It is important to keep in mind that the choice 
of SFR tracers and spatial scales means that different 
studies effectively sample different time scales, thus the 
SF history of any particular galaxy potentially plays an 
important role in determining the result of the measure- 
ment. It is also possible that these differences correspond 
to a spectrum of physical SF mechanisms present in a 
wide range of environments: in that case, the local SF 
law would not be universal. It is, therefore, vital to un- 
derstand the impact of systematics on the measurement 
of the parameters of the local SF law. Whether the lo- 
cal SF law is linear or non-linear has implications for the 
dominant SF mechanisms as well as for modeling efforts. 

The objective of this paper is to explore the molecular 
gas SF law in the galaxy NGC 4254 (M 99), at -500 pc 
and — 1 kpc scales using different SFR tracers. This 
galaxy has been the subject of a number previous SF 
law studies (Kennicutt 1989; Boissier et al. 2003; Ko- 
mugi et al. 2005; Wilson et al. 2009), although not at 
such high resolution. This is a pilot study using observa- 
tions obtained by the Survey Towards Nearby Infrared- 
bright Galaxies (STING; Bolatto et al. in prep.), which 
employs the Combined Array for Research in Millimeter 
Astronomy (CARMA) interferometer. 

In this paper we investigate the impact of various 
methodological aspects related to local SF law study. 
Our study focuses on 1) the use of different SFR tracers 
and the scatter associated with those tracers, 2) the role 
of the diffuse emission, a component of the integrated 
disk emission which is not necessarily related to the star- 
forming regions, and 3) the role of fitting methodologies 
and data sampling strategies in determining the func- 
tional form of the SF law. We do not explore the role 
of variations in the stellar initial mass function (IMF), 
the CO-to-H2 conversion factor, the extinction correc- 
tion, and various other assumptions pertinent to local 
SF law studies. These issues will be addressed in the 
future using other sample galaxies of the STING survey. 

The organization of the paper is as follows. In S}2] we 
present the multi-wavelength data set, including a brief 
description of the data products. A discussion of the 
sky background and the extended diffuse emission (DE) 
is presented in Section SQ] contains our main results 
and general discussions. We compare our results with 
the most recent studies of local SF law in $5] The main 
findings of our study are summarized in fj6] A brief in- 
troduction to NGC 4254, details on the construction of 
various data products such as surface density maps, the 
discussion on the treatment of DE, and the details of 
the sampling and the regression analysis can be found in 
several sections of the Appendix. 

2. DATA 

The target of this study, NGC 4254, is an almost face- 
on (i — 30°) SA(s)c spiral located at a kinematic distance 
of 16.6 Mpc (Prescott et al. 2007). It lies in the outskirts 
of Virgo cluster, — 1.2 Mpc to the north-west from the 
cluster center in projected distance. For our adopted 
distance to NGC 4254, 1" corresponds to 80 pc. The 
optical radius of this galaxy is R25 — 12.1 kpc. See 
appendix \X\ for more information on this galaxy. 

We construct the molecular gas density maps using 
a combination of CO J = 1 — emission obtained by 
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CARMA interferometer for the STING0 survey, and 
single-dish CO J = 2 — 1 observations obtained by the 
Institut de Radio Astronomie Millimetrique (IRAM) 
30 m telescope at Pico Veleta, Spain. These data are 
part of the HERA CO Line Extragalactic Survey (HER- 
ACLES) and were observed and reduced in the same 
manner as the first part of the survey described in Leroy 
et al. (2009). The full survey will be presented by (Leroy 
et al. in prep.). To construct SFR maps of NGC 4254 
we use ultraviolet (UV) images from GALEX Nearby 
Galaxy Survey (NGS; Gil de Paz et al. 2007), and Ha and 
24 fim images from the Spitzer Infrared Nearby Galaxies 
Survey (SINGS3 Kennicutt et al. 2003) archive. 

Here we describe the multi-wavelength data. The ba- 
sic information for the data set is provided in Table [TJ 
For various other properties of NGC 4254 the reader is 
referred to Table 1 of Kantharia et al. (2008). 

2.1. CARMA STING Data 

The interferometric CO J = 1 — map of NGC 4254 
was obtained as part of the STING survey using CARMA 
interferometer. The STING sample is composed of 
northern (S > —20°), moderately inclined (i < 75°) 
galaxies within 45 Mpc culled from the IRAS Revised 
Bright Galaxy Survey (RBGS; Sanders et al. 2003). 
These galaxies have been carefully selected to have uni- 
form coverage in stellar mass, SF activities, and morpho- 
logical types. The survey is complementary to BIMA 
SONG (Heifer et al. 2003) but targeted to have better 
disk coverage, sensitivity and resolution. The details of 
the CARMA STING survey will be published in a forth- 
coming paper (Bolatto et al., in preparation). 

Observations with the CARMA interferometer were 
conducted in the D array configuration during June 2008 
for a total of 8.5 hours. Passband and phase calibration 
were performed using 3C273, and 3C274 was used as a 
secondary phase calibrator to assess the quality of the 
phase transfer and coherence. The absolute flux scale 
for the interferometer was determined by observing Mars. 
At A = 2.6 mm the 6 m (10 m) diameter CARMA an- 
tennas have a primary beam FWHM of 90" (54")which 
defines their effective field of view. The observations 
were obtained using a Nyquist-sampled 19-pointing mo- 
saic pattern that provides an effective field of view of 
100" in diameter. The synthesized beam produced using 
natural weighting has 4.3" FWHM, which is the resolu- 
tion of the resulting map. 

In our frequency setup, the receivers were tuned to the 
CO J = 1 — transition at a rest frequency of 115.2712 
GHz (A = 2.6 mm) in the upper sideband, with the CO- 
BRA correlator configured to have three spectral win- 
dows in each sideband, one in the low resolution 500 
MHz (16 channels) mode, and two partially overlapping 
in the 62 MHz (63 channels) mode. These were placed 
side by side in frequency with an overlap of 10 MHz 
(6 channels). The resulting velocity range covered was 
290 km s -1 with a intrinsic velocity resolution of 2.6 km 
s . The maps were produced with 10 km sec -1 veloc- 
ity resolution. The sensitivity of the interferometric map 
before combination with the single-dish data is a ks 22 

9 http://www.astro.umd.edu/~bolatto/STING/ 

10 http:/ /ssc. spitzer. caltech.edu/legacy/singshistory.html 



mjy Beam 1 in 10 km sec 1 wide channels (see Fig. 

The high angular resolution measurements obtained 
with interferometers filter out the large spatial scales of 
the source, giving rise to the "missing flux" problem. 
As a result, surface densities for extended sources de- 
rived from interferometric data alone may be underes- 
timated. To overcome this limitation it is necessary to 
merge single-dish observation of the object with the in- 
terferometric data. 

The single dish IRAM CO J = 2 - 1 observations at 
A = 1.3 mm of NGC 4254 have ~ 12.5" beam and extend 
beyond the region mapped by STING. The sensitivity of 
the single dish map is a « 27 mK in 2.6 km sec -1 wide 
channels. We have used a gain factor ofT/S ~ 0.14K/Jy 
appropriate for the IRAM 30m single dish telescope and 
converted the units of the CO J = 2 — 1 data cube from 
Kelvin to Jy Beam -1 . 

Comparison of the enclosed fluxes shows that (assum- 
ing thermalized optically-thick CO emission, see below) 
the STING map recovers most of the single-dish flux in 
its inner 60", progressively losing flux beyond that point 
(see Fig. Q]). We converted the CO J = 2 — 1 single dish 
cube to the equivalent CO J = 1 — flux by applying 
a multiplicative factor of 4, following the assumption 
that the (peak) brightness temperature (T m b) is approxi- 
mately the same for the CO J = 1 — and CO J = 2 — 1 
transitions. Before combination the CARMA cube was 
de-convolved using the implementation of the CLEAN 
algorithm in the MIRIAD task "mossdi". We have re- 
binned the velocity channels in the IRAM data cube 
to attain that of the CARMA cube. The combination 
between the CARMA and the IRAM cube was accom- 
plished in the image plane (e.g., Stanimirovic et al. 1999), 
using the MIRIAD task "immerge" . The spatial resolu- 
tion of the combined map is the same as the CARMA 
CO J — 1 — map. We have implicitly assumed a uni- 
form filling factor while combining the two CO maps. 

Validation for the assumption of optically-thick ther- 
malized CO emission comes from the comparison of 
single-dish fluxes. With a beam size of 45", the 
global CO flux for NGC 4254 estimated in the FCRAO 
CO J = 1 - survey is (3000±850) Jy km sec -1 (Young 
et al. 1995). Our estimated flux using the IRAM 
CO J = 2 — 1 observations is ~ 2760 Jy km sec -1 out to 
R25, showing that IRAM reproduces the FCRAO mea- 
surement within ~ 10%. This similarity shows that the 
assumption of identical brightness temperature for the 
1 — and 2 — 1 transitions is reasonable. 

Although the details of the study we present in the 
following sections depend somewhat on the combination 
of the interferometer with the single dish data, the major 
results are very much independent. 

2.2. UV, Optical and Mid-infrared Data 

NGC 4254 has been observed in the near-ultraviolet 
(NUV, 1750-2750 A) by GALEX. To convert the map 
from NUV to FUV we use a morphology dependent color 
correction (FUV - NUV) « 0.46, following Gil de Paz 
et al. (2007). The morphology parameter T = 5 of 
NGC 4254 is obtained from the Third Reference Cat- 
alog (RC3; de Vaucouleurs et al. 1991). To correct 
the FUV map for line-of-sight Galactic extinction we use 
Apuv ~ 8.24 E(B - V) (Weyder et al. 2007), where the 
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Table 1 

Basic Information of the Data Set 



Telescope 


Wavelength 


Pixel 


FWHM 


Sensitivity (ler) 


Sensitivity Unit 


GALEX 


0.2271 nm 


1.5 


5.6 


3.4 ± 1.00 


10~ 15 erg sec -1 cm~ 2 


KPNO 


0.6563 ^im 


0.3 


1.5 


(6.7 ± 1.5) x 10~ 2 


10 — 15 erg sec -1 cm~ 2 


Spitzcr 


24 fim 


1.5 


G.O 


3.3 ± 0.82 


10 — 15 erg sec -1 cm -2 


CARMA 


2.6 mm 


1.0 


4.3 


0.22 


Jy Beam -1 km sec -1 


IRAM 


1.3 mm 


2.0 


12.5 


1.00 


Jy Beam -1 km sec -1 


^SFR 




3.0 


6.0 


0.10 


M Gyr- 1 pc" 2 


Sh 2 




3.0 


0.0 


3.70 


M pc~ 2 



Note. — The pixel resolution and the FWHM of FUV, Ha, and MIR 24 /im maps are in unit of arcsec. The limiting 
sensitivities of the CARMA CO J — 1 — and the IRAM CO J = 2 — 1 observations have been estimated from a single velocity 
channel map and 4 consecutive velocity channel maps, respectively, using a — [a c h\/N] Av, where cr c h is the rms noise in a 
velocity channel map, N is the number of channel, and Av is the velocity resolution. The velocity resolutions for the CARMA 
and the IRAM observations are Av = 10 and 2.6 km sec -1 , respectively. In this table only the surface densities are inclination 
corrected. 



Galactic reddening E(B — V) « 0.04 is estimated from 
Schlegel et al. (1998). The FUV map is converted to 
AB magnitudes using the following formula (Gil de Paz, 
private communication), 



niAB = —2.5 log [counts/sec] + 18.82 



(2) 



The SINGS project public data archive provides cal- 
ibrated and stellar continuum subtracted Ha image of 
this galaxy. Comparing the R-band and Ha images we 
identify foreground stars which are then masked, partic- 
ularly those within the optical diameter. The resulting 
image is then corrected for [N II] AA6548, 6583 forbidden 
line emission and the transmission curve of the Ha filter, 
using the factors obtained for this galaxy by Prescott et 
al. (2007). 

The SINGS archive also provides the mid-infrared 
(MIR) 24 /im image, which s a scan map taken with the 
MIPS instrument on board the Spitzer Space Telescope 
(Rieke et al. 2004). The MIPS data were processed using 
the MIPS Instrument Team Data Analysis Tool (Gor- 
don et al. 2005). No stellar masking was necessary for 
the MIR map of NGC 4254. Basic information of these 
images are given in Table [T] 

2.3. Data Products 

The spatial resolution of our study is limited by the 
point spread function (PSF) of the MIR data, which has 
a FWHM of 6". This angular scale corresponds to a 
physical length of ~ 480 pc in the disk of NGC 4254. We 
should note that, although we approximate it as a Gaus- 
sian, the mid-infrared PSF is complex. It has prominent 
first and second Airy rings, with the second ring stretch- 
ing out to ~ 20". Nevertheless, approximately 85% of 
the total source flux is contained within the central peak 
with FWHM of 6" (Engelbracht et al. 2007). 

The higher resolution Ha and CO images were 
Gaussian-convolved to have the same resolution and sam- 
pling as the MIR image. In both cases, the convolution 
and regridding used the AIPS packageEJ No convolution 
was necessary for the FUV image, since it has a resolu- 
tion similar to the MIR (see Table 1). For our high- 
resolution analysis we regrid the images to 3" pixels to 
Nyquist-sample the PSF. 

11 The Astronomical Image Processing System (AIPS) has been 
developed by the National Radio Astronomy Observatory (NRAO) 



Figure [2] shows the FUV, Ha, and 24 /jm images of 
NGC 4254 used to construct the SFR maps in logarith- 
mic color scale. The black contours correspond to the 3 
and 12 Jy Beam -1 levels of the CO J — 1 — map. Fig- 
ure [2] shows the striking similarity between the distribu- 
tions of hot dust (T ~ 100 K) traced by the 24 /im emis- 
sion and the cold molecular gas traced by CO J = 1 — 
in NGC 4254. Interestingly, the CO J = 1 - map 
shows highly symmetric north and south spiral arms, 
whereas the northern spiral arm is not prominent at op- 
tical and FUV wavelengths. 

Construction of various data products such as molec- 
ular gas and SFR surface density maps and the asso- 
ciated error maps are summarized in Appendix [B] and 
ICl The surface density maps are de-projected using the 
inclination angel i ~ 30°. In this study we construct 
four different SFR surface density maps from combin- 
ing FUV and 24 /im (Spuv+24/im), extinction corrected 
Ha (Sh q ), observed Ha and 24 /im (SHa+24/xm), and 24 
/im (£24 /im) following various prescriptions in the litera- 
ture. All four of SFR surface density maps are expressed 
in units of M Q Gyr -1 pc -2 . The limiting la sensitiv- 
ity in surface density varies from map to map where the 
Ha + 24 /im map has the highest rms sensitivity (~ 0.1 
M Q Gyr -1 pc -2 ) among all four SFR tracer maps. We 
adopt this value as the limiting sensitivity for all Esfr 
maps. 

The inclination corrected limiting surface densities cor- 
responding to the sensitivity limits (Table [1]) of the 
CARMA interferometer and the IRAM single dish maps 
are ~ 5.3 M Q pc -2 and ~ 0.6 M Q pc -2 , respectively. As 
mentioned in section [2~T| we apply a multiplicative factor 
of 4 to the sensitivity limit of the IRAM single dish map 
to derive the limiting surface density. We combine the 
interferometric and single-dish maps to create the molec- 
ular gas surface density map (£h 2 )- This combined map 
finally convolved with a Gaussian beam to the obtain the 
spatial resolution of 6". The typical lcr sensitivity of this 
Sh 2 map is ^3.7 M Q pc -2 (inclination corrected). We 
study the molecular gas-SFR surface density relation for 
Sh 2 ~ 10- 245 M pc -2 . 

2.4. Sampling and fitting considerations 

Since one of the goals of this study is to explore how the 
functional form of the SF law depends on the treatment 
of the data, we analyze the images using 1) pixel anal- 
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ysis, incorporating all the data above a signal-to-noise 
cut, 2) aperture analysis, where we average over circu- 
lar apertures selecting bright regions, and 3) azimuthally 
averaged annuli, with a width of 500 pc. 

In local SF law studies, especially for normal star- 
forming galaxies such as NGC 4254, the dynamic range 
probed by the molecular gas is rather small. The ob- 
served dispersion in the SFR tracer, on the other hand, is 
usually quite large depending on the selection of the SFR 
tracer. Due to this characteristic of the gas-SFR surface 
density relation, the determination of the functional form 
of this relation depends critically on statistical method- 
ologies and fitting procedures. We describe the sampling 
and fitting strategies in detail in the appendix sections 
O and El 

3. DIFFUSE EXTENDED EMISSION 

Several components contribute to the total emission 
in galaxies. Images contain emission from backgrounds 
or foregrounds, which are not physically related to the 
galaxy. Besides the emission of the localized SF, they 
also contain diffuse components that are extended over 
the entire disk and not necessarily associated with SF 
activity, which we discuss in the following section. Since 
the calibration of SFR tracers is frequently performed in 
star-forming regions, it may be important to remove the 
contribution from diffuse components to the brightness 
distribution before interpreting it in terms of a SFR. 

The CO distribution of a galaxy can also contain dif- 
fuse emission not necessarily associated with the indi- 
vidual star-forming regions. A collection of unresolved 
small molecular clouds, in particular Taurus-like clouds 
in the Milky Way with masses M - 10 4 M , will fall 
below our detection threshold as individual entities but 
would contribute to the diffuse extended emission. It 
should be removed if those clouds do not host massive 
SF contributing to the SFR tracers. 

3.1. Diffuse Emission in the SFR Tracers 

The DE is ubiquitous in the UV, Ha and 24 /jm maps 
and it spreads out across the disk over a few tens to 
hundreds of parsecs from the clustered OB association 
and resolved H II regions. The origin of this emission an 
active area of research over the past four decades (see 
Monnet 1971; Reynolds et al. 1971; Haffner et al. 2009 
and references therein for diffuse emission from H II re- 
gions; see Witt 1968; Hayakawa & Yoshioka 1969; Meurer 
et al. 1995; Pellerin et al. 2007 and references therein for 
diffuse UV emission; see Dale et al. 2007; Draine et al. 
2007; Verley et al. 2007, 2009 for diffuse 24 fim emission). 

The extended Ha emission originates in diffuse (n e ~ 
0.1 enr 3 ) warm (T-8000 K) ionized gas (DIG), which 
is analogous to the warm ionized medium (WIM) in our 
own galaxy. This gas has a volume filling factor of —0.25 
and scale height of — 1 kpc and is a major component 
of the interstellar medium of the Galaxy (see Reynolds 
1991, 1993 for reviews). In the Milky Way the DIG con- 
tributes ~ 10 — 15% to the total H II emission. For ex- 
ternal galaxies, however, observational evidences suggest 
that the DIG may contribute a substantial fraction (~30- 
60%) to the total emission, fairly independent of galaxy 
Hubble type, inclination, and SF properties (Hunter & 
Gallagher 1990, 1992; Rand et al. 1990; Walterbos & 
Braun 1994; Kennicutt et al. 1995; Hoopes et al. 1996, 



2001; Ferguson et al. 1996; Hoopes & Walterbos 2000; 
Wang et al. 1997; Greenawalt et al. 1998). Using spec- 
tral data Blanc et al. (2009) very recently reported 11% 
DIG contribution in the central 4 kpc in M51. They 
also find that the DIG makes 100% of the total emission 
coming from the interarm regions of this galaxy. 

Is the emission from the DIG a tracer of SF ? On the 
one hand, the large energy output and the morphologi- 
cal association of the DIG with H II regions have been 
used to argue that early type OB stars in H II regions are 
the sources of this diffuse emission. Leakage of ionizing 
photons from porous H II regions has been invoked to 
explain the widespread distribution of this component. 
This requires the interstellar medium (ISM) to have low 
FUV extinction along certain lines-of-sight, allowing a 
large mean free path for these photons likely through in- 
terconnecting ionized bubbles (Tielens 2005, Scon 2009). 
On the other hand, observational studies suggest that the 
DIG may not be entirely associated with the early-type 
massive OB stellar clusters in H II regions. A popula- 
tion of late-type field OB stars (Patel & Wilson 1995a, b; 
Hoopes et al. 2001) or supernovae shocks may also pro- 
vide energy to the DIG (Collins & Rand 2001; Rand et al. 
2007). The relative contribution of each of these sources 
to the DIG energy balance is uncertain. Recent numer- 
ical simulations, however, suggest that the contribution 
from H II regions to the energy budget of the DIG could 
be -30% or less (Seon 2009). 

Because mostly non-ionizing photons contribute to it, 
the diffuse UV continuum emission has an origin different 
from that of the DIG. In starburst galaxies, Meurer et 
al. (1995) found that about 80% of the UV flux at 2200 
A is produced outside clustered OB associations and it 
has an extended character. Popescu et al. (2005) sug- 
gested UV light scattered by dust as the possible origin 
of the diffuse UV emission. Tremonti et al. (2001) and 
Chander et al. (2003, 2005), however, have noted that 
for starburst galaxies the spectral UV lines from clusters 
are different from those in the inter-cluster environment. 
Their studies show that the UV stellar signature in clus- 
ters is dominated by O-type stars, while the inter-cluster 
environment is dominated by less massive B-type stars. 
Late type OB field stars were also suggested by Hoopes 
et al. (2001) as the origin of the diffuse UV emission in 
normal spirals. In a recent study, Pellerin et al. (2007) 
find that -75-90% of the UV flux is produced by B-type 
field stars in the disk of the barred spiral NGC 1313. 
These studies suggest that B-type field stars are the ma- 
jor source of non-ionizing UV emission in galaxies, with 
a much smaller contribution from scattered OB cluster 
light. This implies that the SF history has an impor- 
tant role in determining the ratio between the diffuse 
UV continuum and that arising in compact OB associ- 
ations. Late B-type stars are longer lived (— 100 Myr) 
and less massive (~ 5 — 2OM ) than O-type stars (age 
—10 Myr, mass > 20 M Q ), with the latter types mostly 
residing in clustered associations. 

The 24 yum continuum emission also has a diffuse com- 
ponent associated with it. In galaxy disks, 24 /im dust 
emission is frequently found near discrete H II regions 
(Helou et al. 2004). This extended 24 fim emission is due 
to small dust grains out of equilibrium with the radia- 
tion field, for which single-photon events produce large 
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Table 2 

Diffuse Fraction at Various Filter Scales 



Number 


Filter Width 


Filter Width 




Diffuse Fraction 






(arcscc) 


(kpc) 


/fuv 


/dig 


/mir 



I 


75 


6.03 


0.72 


0.55 


0.68 


If 


90 


7.23 


0.69 


0.50 


0.62 


III 


105 


8.44 


0.66 


0.45 


0.56 


IV 


120 


9.64 


0.62 


0.40 


0.50 


V 


135 


10.85 


0.59 


0.35 


0.44 


VI 


150 


12.06 


0.55 


0.30 


0.38 


VII 


165 


13.27 


0.52 


0.25 


0.34 


VIII 


180 


14.47 


0.48 


0.21 


0.30 


IX 


195 


15.68 


0.45 


0.18 


0.25 


X 


210 


16.88 


0.42 


0.15 


0.22 


XI 


225 


18.10 


0.39 


0.11 


0.19 



Diffuse fractions as a function of galacto-centric radius that area obtained from different SFR maps arc shown in Fig. I14l in the appendix 



Note. — 

m 

temperature excursions (Desert et al. 1990). In addition 
to this localized emission, 24 /jm sources are surrounded 
by a diffuse component associated with overall galaxy 
profile and internal structure such as spiral arms (Helou 
et al. 2004; Presscott et al. 2007; Verley et al. 2007, 
2009). The old stellar population is thought to be respon- 
sible for such component, which comprises ~ 30 — 40% 
of the total thermal dust emission in the central regions 
to ~ 60 — 80% of the integrated emission in the extended 
disk (Draine et al.2007; Verley et al. 2009). 

Understanding the nature and sources of DE is of great 
importance in studies of SF. Kuno et al. (1995) and Fer- 
guson et al. (1996) discussed the role of the diffuse com- 
ponent when deriving the SFR based on Ha emission. An 
assessment of the magnitude of DE contribution is neces- 
sary in order to use the SFR tracers derived from FUV, 
Ha, and MIR 24 /im dust emission. Thus, it is interest- 
ing to study the contribution of the DE in these tracers 
and its consequences on spatially resolved molecular gas 
SF law studies. The contribution from DE is most im- 
portant in the low surface brightness regime, where it can 
flatten the power law index of the SF law if unaccounted 
for. 

3.2. Diffuse Emission in the Molecular Gas 

In the Milky Way most CO emission is associated with 
GMCs, which have a top heavy mass function (most of 
the mass is in the most massive GMCs; Solomon et al. 
1987). Similar "top heavy" CMC mass functions are 
observed in most Local Group galaxies with the excep- 
tion of M33 (Engargiola et al. 2003). Even in the Milky 
Way, however, there exists a population of CO-emitting 
molecular clouds that are considerably more diffuse and 
have lower masses and column densities than GMCs that 
host massive SF. Examples are the high latitude clouds 
(Magnani et al. 1985), with typical column densities of 
Nh2 ~ 10 21 cm~ 2 and very low SF activity. 

We do not know whether the GMC mass function 
in NGC 4254 is top heavy or bottom heavy. Even if 
it is top heavy, at 500 pc resolution our observations 
(la S H2 ~ 3.7 M pc~ 2 ) will be sensitive only to GMCs 
with masses M > 10 6 M Q as distinct entities. Local- 
ized, lower mass GMCs will be blended together and will 
appear as a blurred diffuse emission background in the 
galaxy disk. If these lower mass GMCs host no massive 
SF contributing to the selected SFR tracer they should 
not be included in the determination of the SF law, oth- 



erwise their inclusion will artificially steepen the power 
law index of the SF law. 

3.3. Removing Diffuse Emission 

We will discuss first the treatment of the DE in the 
SFR tracers. The procedure for removing the DE in the 
molecular gas map is very similar, and is discussed in 

El 

There is no standard procedure to separate integrated 
flux into emission from discrete star-forming regions and 
an extended diffuse component. Various criteria have 
been adopted in previous studies to distinguish between 
the discrete H II regions and the DIG in the nebular 
emission map. For example, forbidden line ratio of 
[S II] AA6716, 6731 to Ha (Walterbos & Bruan 1994), 
the equivalent width of Ha line (Veilleux et al. 1995), 
surface brightness cut (Ferguson et al. 1996), unsharp 
masking (Hoopes et al. 1996), H II region luminosity 
function (Thilkcr et al. 2002). Although the forbidden 
line ratios are powerful probes to separate the compo- 
nents unambiguously, in absence of such information two- 
dimensional (2D) spatial filtering techniques such as un- 
sharp masking is the most robust method to obtain in- 
formation about DE. 

We use unsharp masking to model and remove the DE, 
taking advantage of its extended nature. Several authors 
have used unsharp masking to separate DE and discrete 
H II regions in Local Group galaxies (Hoopes et al. 1996; 
Greenawalt et al. 1998; Thilker et al. 2005). Our ap- 
proach, slightly different than these studies, is simple, 
easy to implement, and avoids an ad hoc surface bright- 
ness cut used in other studies. The process involves cre- 
ating a smoothed or blurred image produced by a 2D 
moving boxcar kernel (middle panels, Fig. [3]) and then 
the subtraction of the smoothed map from the original 
image (see Appendix [Fl for a more detailed explanation). 
It is possible to use other kernels (Gaussian, Hanning, 
etc), but boxcar is the computationally simplest. The 
resulting final image (lower panels, Fig. [3]) has a reduced 
contribution from the background as well as the DE, 
since most of it is contained in the spatially smoothed 
map. Indeed deep Ha images of the distribution of DIG 
in the local group members, such as M 33 and M 81, 
show it to be quite smooth (Greenawalt et al. 1998). 

Since we are using multi-wavelength SFR tracers, it is 
important to understand the nature and brightness of the 
DE as observational studies suggest that these properties 
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depend on the wavelength studied. For example, studies 
of the fraction of DE at 24 fan in the disk of M 33 find 
/mir. ~ 60 — 80%, increasing radially outward (Verley 
et al. 2009). By contrast, the DIG shows the opposite 
trend, with /dig ~ 60% at the center and decreasing to 
almost zero towards the outer disk. It is generally near 
40% across the disk (Thilker et al 2005). The diffuse 
fraction in FUV shows a remarkably flat profile /fuv ~ 
0.65 (Thilker et al. 2005). 

The crucial aspect of unsharp masking is the choice of 
the size of the median filter kernel. The filtering kernel 
size affects the fraction of the total emission of the orig- 
inal map contained in the smoothed image, /de, with 
the larger fractions in the smooth or diffuse component 
corresponding to the smaller filtering kernel sizes. In the 
subsequent analysis we will use /de to refer the diffuse 
fraction in general. To refer to the diffuse fraction in the 
FUV, Ha, and 24 /jm we use the notation /fuv, /dig, 
and /mir, respectively. 

The details of the unsharp masking process can be 
found in the appendix [F] In NGC 4254 we explore a 
number of filter sizes in each SFR tracer, carrying out 
our analysis for each case (see Table [2]). The diffuse frac- 
tion as a function of filter scale is shown in the panel (D) 
in Fig. [TJ] At a given filter scale the diffuse fraction 
is different for different SFR tracers. For convenience of 
presentation, therefore, we use /dig as the reference DE 
since it is widely known in the astronomical community. 
We will refer the DE as the dominant, significant, sub- 
dominant, and negligible part of the total disk emission 
in the Ha map for /dig > 50%, -30-50%, -10-30%, and 
< 10%, respectively. 

It should be noted here that the noise in the spatially 
smoothed map is negligible compared to the original im- 
age. The subtraction of the smoothed map, therefore 
does not change the noise properties of the original map. 

4. RESULTS AND DISCUSSION 

We begin our analysis presenting azimuthally averaged 
radial distributions of surface densities to demonstrate 
their spatial variations. Fig. [4] shows the derived surface 
densities of the SFR tracers as well as CO gas maps at 
/de = 0. Each panel shows the la dispersion in the 
radial distributions out to O.6R25, where R25 — 12.5 kpc 
is the optical radius of NGC 4254. 

The figure also shows the radial profiles of molecular 
gas obtained from the interferometer-only observations 
and from the combined single-dish/interferometer data. 
Both maps are in excellent agreement within the central 
80" (~ 6.5 kpc), showing that the interferometric data by 
itself accurately traces the distribution of cold molecular 
gas in the inner — 6.5 kpc of this galaxy. The combined 
map, on the other hand, traces better the low surface 
brightness CO in the outer regions of the disk 

4.1. The SF Law using Pixels and Apertures 

In this section we present our results for the case when 
DE is subtracted only from the SFR tracer maps. A 
more general case which addresses DE in both the SFR 
and CO maps will be presented in the following section. 
The Nyquist sampling rate at the fixed 2.5cr cut results 
in — 800 approximately independent pixels for both gas 
and SFR surface density maps at the dominant diffuse 
fractions. The number of pixels increases to — 950 when 



the contribution of DE is sub-dominant or negligible, be- 
cause fewer pixels fall below the signal-to-noise cut after 
DE subtraction. 

We show our results for the Esfr - ^h 2 relation in Fig. 
[6l The figure shows the gas-SFR surface density relation 
for pixel sampling at various /de- The gray scale rep- 
resent the two dimensional histogram of the frequency 
of points, and the contours are placed at 90%, 75%, 
50%, and 25% of the maximum frequency. The diago- 
nal dotted lines represent lines of constant SF efficiency 
(e), or constant molecular gas exhaustion timescale (rdep) 
with values of e = 1%, 10%, and 100% corresponding to 
exhaustion times of Td ep = 10, 1, and 0.1 Gyr respec- 
tively. The filled circle in each panel represents the disk 
averaged surface densities measured within R25 before 
unsharp masking. Within the range of diffuse fractions 
probed the £24 /im - £h 2 relation shows the tightest cor- 
relation whereas the Ehq - £h 2 relation shows the largest 
scatter. We compute the linear Pearson correlation coef- 
ficient (r) for these two relations in the range of explored 
diffuse fractions, finding r — 0.55 — 0.7 for the former and 
r — 0.4 — 0.55 for the latter. The observed dispersion is 
<7i — 0.3 dex for £24 /im - £h 2 an d ~ 0.5 dex for Eh q - 
^h 2 - 

Note that the scatter in the SF law is substantially 
lower when no DE is subtracted from the total emission 
of the SFR tracers. Furthermore, since the DE is propor- 
tionally more important in fainter regions, its subtraction 
increases scatter in the gas-SFR relation mostly at low 
surface densities. For the same reason, removing the DE 
steepens the SF law. 

The results of aperture sampling are shown in Fig. [7] 
11 

for a 105 unsharp masking kernel. The distributions of 
points are overlaid on the contours obtained from the 
pixel analysis. By construction the apertures sample 
mostly the high density regions, and the overall agree- 
ment in these regions is excellent between the pixel and 
the aperture analysis. The lack of the low surface bright- 
ness tail along the vertical axes, however, has important 
implications for the slope of the SF law, as we discuss 
next. 

We show the measurements from different bivariate re- 
gression methods in Fig. [8] The figure highlights the 
£Ha+24^m " ^ H 2 relation for pixel sampling at 105" fil- 
ter scale. The scale corresponds to a case of significant 
to dominant diffuse fraction, /dig=0-45 and /mir=0.56. 
The figure shows that the FITEXY method yields the 
shallowest slope. For this gas-SFR surface density rela- 
tion the power-law index is in the range N mo i — 1.3 — 1.6 
for pixel sampling, with intrinsic scatter o\ — 0.3 — 0.4 
obtained from the fits. All fitting methods yields system- 
atically shallower slope (N mo i — 0.8 — 1.2) and smaller 
scatter (<r; — 0.1 — 0.3 for both aperture sizes) for aper- 
ture sampling. 

Our exploration of the range of results obtained for 
the SF law by using different methodologies is summa- 
rized in Fig. |9] The figure shows the dependence of the 
power-law index of the SF law on the subtracted diffuse 
fractions. The vertical bars associated with each point 
are purely methodological, and show the range of slopes 
obtained from applying the different fitting methods de- 
scribed in appendix [E] We refer to them as the method- 
ological scatter, <7 mo i,fit- The mean (Nmo^gt) obtained 
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from averaging the results of the three fits for each /de 
is shown with filled circles. The panels in this figure il- 
lustrate how the measurement depends not only on the 
chosen SFR tracer, but also on the type of analysis and 
on the treatment of the DE. The linearity of the func- 
tional form of the molecular SF, in particular, depends 
on the amount of DE assigned to either axis. It is impor- 
tant to notice, however, that this slope change is driven 
by the lowest surface brightness regions of the disk. In 
the high surface-brightness regions sampled in the aper- 
ture analysis the choice of /de is unimportant, and a 
unique slope is consistent with the data for any (rea- 
sonable) amount of DE. In these regions the SF law is 
approximately linear, although its precise value depends 
on the SFR tracer. 

We observe a direct relation between the slope of the 
SF law and the magnitude of the DE subtracted in 
the pixel analysis (left panels, Fig. [9J. For a domi- 
nant diffuse fraction (/de ^ 50 — 60% of the total disk 
emission), all the resolved SF law relations in the pixel 
analysis show systematically the steepest power-law in- 
dices, Nmoi^t ~ 1-4 — f.7. For sub-dominant to neg- 
ligible diffuse fraction (/de ^ 30% of the total disk 
emission), however, the slope clearly becomes shallower, 
NmoKfit ~ 1 — f.2. Thus higher /de corresponds to a 
steeper power-law index. This is only observed in the 
pixel analysis, which contains the low surface brightness 
regions. Furthermore, the scatter in the results yielded 
by the different fitting algorithms is also a monotonically 
increasing function of the amount of DE subtracted. This 
methodological scatter is driven by the corresponding in- 
crease in the scatter of the low surface brightness pixels, 
which have a very broad distribution for large /de- 

For the aperture analysis the fitted power indices are 
systematically shallower than for pixel analysis, and ro- 
bust to the choice of /de and the aperture size (right 
panels in Fig. [9]). For a dominant diffuse fraction 
(/de }t 50%), the measurements of the SF law slope 
are consistent with Nmo^fit ~ I. For small diffuse frac- 
tions (/de S 20%) diffuse fraction, the fitted power- 
index is in the range Nmoi^t ~ 0.8 — f . Furthermore, the 
Ha corrected for azimuthally averaged extinction tends 
to consistently have the steeper slopes (and the highest 
methodological scatter). Although with this data sam- 
pling the slope still flattens monotonically with the re- 
duction of the amount of subtracted DE,_the dependence 
on /de is very weak and the variation in Nmo^gt is within 
the scatter of the different fitting methods. 

The power law index is slightly steeper for the larger, 
1 kpc diameter apertures at the dominant and signifi- 
cant diffuse fractions where the fitting procedures diverge 
more from one another. This is likely due to a combi- 
nation of the fact that the larger apertures encompass 
some area of low surface density material, and to the re- 
duction in the number of data points by a factor of ~ 3. 
With fewer data the regression analysis becomes highly 
sensitive to the distribution of points. For sub-dominant 
to negligible diffuse fractions (/de ^ 0.3), however, the 
power law index of the local SF agrees well for both aper- 
ture sizes. 

Tables |4] - [6] show the derived parameters for median 
filter of size 75", 105" and 180". The row represented 
by dash in each table show the parameter for /de=0, 



i.e., when the filter size is the same of the entire map. 
The filter widths are chosen to show the representative 
cases of the dominant, significant, sub-dominant, and 
negligible diffuse fractions. At a given filter scale the 
estimates from the OLS bisector, FITEXY, and LINEX- 
ERR methods are shown by the top, middle, and, bottom 
row, respectively. The quoted error in each parameter 
comes from bootstrap sampling of 1000 realizations of 
data points. 

4.1.1. The SF Law in Annuli 

Many of the early resolved studies of the relation be- 
tween gas and star formation in galaxies analyzed the 
data using azimuthal averages (e.g., Wong & Blitz 2002). 
Following the procedure similar to that discussed in ap- 
pendix[D]to select the common regions from the £h 2 and 
Ssfr maps, we also explore the SF law for azimuthally- 
averaged radial profiles (Fig. ITO)) . Sampled in this man- 
ner, the_ functional form of the SF law in NGC 4254 is 
linear (Nmo^fit ~ 1) for /de = 0, and approximately lin- 
ear (N m oi,fit ~ 1 — 1.2) in the range of diffuse fractions 
studied. The linear form stems from the fact that the 
azimuthal averages are dominated by the high surface 
brightness regions, and there is no "extended tail" of low 
SgFR points steepening the fit to the distribution. Since 
the data have low dispersion all fitting methods yield 
consistent results. Table [7] shows the fitted parameters 
derived from the OLS bisector method in all SFR tracers 
at the diffuse fractions highlighted in Fig. [TUl 

4.2. Dispersion in the Relations 

The intrinsic dispersion (ai; see appendixlEl) in the gas- 
SFR surface density relations at various diffuse fractions 
is shown in Fig. 1111 The figure shows the mean (<7i,fit) 
and the scatter ("dispersion of dispersion") obtained by 
the three regression methods. The SFR obtained from 
24 /jm displays the tightest correlation with the molec- 
ular gas, among all tracers (<7; ~ 0.1 — 0.3 dex). This is 
likely due to a combination of two effects: 1) young star- 
forming regions that are still embedded in their parent 
clouds will emit brightly at 24 fim. Bright Ha emission 
will only happen when the HII is older and the parent 
cloud is at least partially cleared (see also Helou et al. 
2004; Relano & Kennicutt 2009). And, 2) by its nature, 
this SFR tracer does not need to be corrected by extinc- 
tion. The spatial correspondence between the 24 /an and 
CO maps is striking (Fig. 2). 

The extinction-corrected Shc, on the other hand, 
shows the largest scatter (ct; ~ 0.3 — 0.6 dex) of all trac- 
ers, with or without unsharp masking. This is because 
the extinction correction is azimuthally averaged, and 
it does a poor job at correcting any position although 
it yields the correct result in a statistical sense. Using 
one galaxy- wide correction factor to remove the [N II] 
AA6548, 6583 forbidden line emission from the Ha map 
is also another potential (likely minor) contributor to the 
scatter, since it may well vary with the position. 

Due to its large dispersion, the results for Ehq from 
the different regression methods differ substantially from 
one another (<7 mo i,fit ~ 0.4). By contrast the com- 
bined Ha + 24 fxm tracer, which applies the same un- 
derlying extinction correction locally, yields a tighter 
correlation (aifit ~ 0.1 — 0.5 dex) and a flatter slope 
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(N mo i, fit ~ 1 - 1.6). The composite FUV + 24 fim yields 
very similar results to Ha + 24 /im. The observed scat- 
ter also becomes somewhat smaller for larger apertures 
(dashed-dot line in the right panels of Fig. [TT|) . partic- 
ularly for Ha which clearly benefits from averaging over 
larger regions. 

4.3. Diffuse CO Emission 

So far we have only considered the effect of diffuse 
emission, possibly unrelated to recent massive SF (on 
the vertical axis of the Ssfr- £h 2 plots). Should we 
be also concerned about analogous effects in the hori- 
zontal axis f £|3.2[) ? To explore the effect of removing a 
diffuse extended component from the CO distribution we 
use smoothing kernels of the same size in both the hor- 
izontal and vertical axes in molecular gas-SFR surface 
density plots. 

Figure [12] highlights the results for FUV + 24 /xm and 
Ha + 24 /im using pairs of lines to illustrate the method- 
ological dispersion. The thin solid and dashed lines show 
the dependence of N mo i on /de when both axes are sub- 
ject to unsharp masking, at two different values of the 
signal-to-noise threshold for including points. To serve 
as comparison, the thick solid lines represents the case 
when only the SFR tracer maps have undergone un- 
sharp masking. In most of our analysis we have only 
included points where the gas surface density map is 
Eh 2 > 2.5ct - 9.2M pc~ 2 (dashed lines in Fig. [Hjl . 
To explore the effects of this threshold on the analysis 
we also plot the results for £h 2 > \o ~ 3.7M Q pc~ 2 
(thin black lines). The threshold for the SFR maps is 
kept at 2.5<r. 

This figure shows that our attempt at removing a dif- 
fuse molecular component in NGC 4254 has only very 
mild impact on the results of the analysis, and only for 
the pixel analysis. The results derived from the aper- 
tures in the high surface brightness regions are essen- 
tially unchanged. Interestingly, the consistency between 
the different fitting methods is better than in the case 
where only the DE in SFR is removed. This is likely be- 
cause errors in /de subtraction smear the data along the 
main relation, rather than only in the vertical direction 
in gas-SFR surface density relation. 

Lowering the threshold after unsharp masking the CO 
produces somewhat flatter slopes at higher /de, while 
increasing the dispersion of the results. For example, the 
slope is approximately unity below /mir ~ 0.4 for a £h 2 
cutoff value of 7.4 M Q pc~ 2 , while for a threshold of 10 
M Q pc -2 it would be unity only below /mir ~ 0.25. The 
power-law index remains unchanged at the extremes of 
/de- The results for other two tracers are qualitatively 
similar to those presented in Fig. [12] 

The dispersion in the gas-SFR surface density relation 
systematically goes down when both variables are subject 
to unsharp masking. We find ~ 10 — 20% reduction in 
the scatter depending on the SFR tracer. The fitting 
methods tend to converge with one another because of 
the reduction in the scatter in the range /de $0.3, which 
is clearly evident in Fig. [12] 

4.4. Goodness of Fit 

How well does the power law functional form represent 
the SFR and molecular gas surface density relationship 



analyzed in this study? A measure of the goodness-of- 
fit of a model is to derive the x 2 statistic based on the 
least square method (Deming 1943). The best fit lines 
provided by the FITEXY estimator have reduced-^ 2 ~ 
1 with probability ~ 0.35 — 0.5. This is achieved by 
iteratively adjusting the error along the Y-axis, cr 2 = 

°m + °fi where o~i is the intrinsic scatter in the gas-SFR 
surface density relation and <r m is the measurement error. 

A graphical alternative to evaluate the goodness-of- 
fit is to test the normality of residuals. The residuals 
are the deviations of observational data from the best fit 
line. We perform this test for the fitted lines produced 
by all three estimators. At large /de the distributions 
of residuals for various SFR tracers are approximately 
consistent with a normal distribution, and they become 
more so when /de decreases. Our analyses suggest that 
the observed relation between the molecular gas and SFR 
surface densities in NGC 4254 is consistent, at least to 
first order, with the power law form. 

4.5. Star Formation Efficiency 

The star formation efficiency (SFE) is a convenient, 
physically motivated way to parametrize the relationship 
between molecular gas and SFR. The SFE has been de- 
fined in various ways in literature. For example, it is 
defined as the ratio of the produced stellar mass to the 
total gas mass. This definition is more commonly seen 
in the Galactic studies (Myers et al. 1986) but also used 
in galaxy modeling (Vazquez- Semadeni et al. 2007). For 
extra-galactic studies, the molecular gas SFE is usually 
defined as (Young & Scoville 1991; McKee & Ostriker 
2007), 

SFE = e =^™=A [Eh/- 1 - 1 . (3) 

^H 2 

The inverse of the SFE is considered as the gas deple- 
tion timescale, Tdep = SFE -1 . This parameter is used to 
discern between the starburst and normal star-forming 
galaxies (Rowand & Young 1999). For starburst galaxies 
the typical depletion time is hundreds of Myr whereas 
normal star-forming galaxies have depletion timescale of 
- 2 Gyr (Bigiel et al. 2008; Leroy et al. 2008). We adopt 
the definition in Eq. [3] for consistency with the studies 
of Bigiel et al. (2008), Leroy et al. (2008), and Blanc et 
al. (2009). 

A linear functional relationship between the SFR and 
molecular gas surface densities implies that the efficiency 
to turn molecular gas into stars is constant across the 
disk. It also implies the gas consumption time is similar 
for both massive and low mass clouds. A linear molec- 
ular gas SF law is consistent with the scenario in which 
GMCs turn their masses into stars at an approximately 
constant rate, irrespective of their environmental param- 
eters. Observations suggest that GMCs properties are 
fairly uniform across galaxies (Blitz et al. 2007; Bolatto 
et al. 2008). A non- linear molecular gas SF law, however, 
implies that gas is turned into stars at a faster rate at 
higher surface densities. 

Figure [13] shows the disk averaged molecular SFE 
(Tdep) as a function of diffuse fraction. We derive the SFE 
finding the ratio of SFR to molecular gas for each pixel 
or aperture in the map, and plot the average with error 
bars computed from the standard deviation. The SFE 
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Table 3 

Disk Averaged Parameters 



Par am. 


DE 


FUV + 24 ^m 


Ha 


Ha ± 24 /Ma 


24 jxra 


Unit 


SFR 


/de=0 


4.9 ±0.1 


2.4 ±0.1 


3.5 ±0.1 


4.1 ±0.2 


M Q yr' 1 


SSFR 


/de=0 


13.3 ± 1.2 


13. 7± 1.2 


11.1 ± 1 


10.5 ± 1 


M Q Gyr" 1 pc~ 2 



r dep l max. / DE 4.3 ±1.2 5.5 ± 1.9 5.4 ± 1.8 4.7 ±1.5 Gyr 

T dep 2 max. /de 2.0 ±0.5 2.5 ± 0.8 2.4 ± 0.7 2.2 ± 0.6 Gyr 

r dep /de=0 1.7 ±0.3 3.3 ±0.8 2.2 ± 0.4 2.0 ± 0.3 Gyr 

Note. — The parameters are estimated within R25 from the maps with zero subtraction of DE, i.e., maps with /de=0. The 
Tdep 1 and Tdep 2 are estimates at maximum /de when 1) only the SFR tracer map, and 2) the SFR tracer and gas maps both 
are subject to unsharp masking. The total molecular gas mass within R 25 is, M H2 ,tot = (5.1 ±0.9) x 10 9 M Q . The disk average 
gas surface density is, £h 2 = (46 ± 3) Mq pc -2 . 



is robustly determined and independent of /de when 
/de J$ 0.3 for all SFR tracers. This is within the range 
of ionized gas /de observed in the Milky Way and Local 
Group galaxies (e.g., Thilker et al. 2002), and also in 
recent spectroscopic determinations in the central region 
of M 51 (Blanc et al. 2009). At higher / DE the SFE 
changes (becomes lower) by up to a factor of 1.5 — 3 at 
the extreme. The global efficiency is approximately inde- 
pendent of /de (within 40%) when both the SFR and the 
gas map have a diffuse component removed (gray points), 
reflecting the fact that the SF law is linear in that case. 
At a given /de, the global SFEs derived from our four 
tracers are approximately consistent with one another, 
although the SFE obtained from Hcv is only marginally 
so (see Table El). The SFE in NGC 4254 is essentially 
independent of radius up to R ~ O.5R25 ~ 6 kpc. Inside 
O.IR25 the SFR derived from the extinction-corrected Hcv 
image does not agree with the other SFR tracers, likely 
because of a problem with the extinction correction. 

Figure [13] also shows that from sub-dominant to negli- 
gible /de the molecular SFE in NGC 4254 is fairly typ- 
ical of large spirals (see Table The disk averaged 
Tdep ~ 2.2 ± 0.4 Gyr estimated from EHa+24/xm map at 
/de = is in good agreement with the Tdep ~ 2.1 Gyr 
measurement by Wilson et al. (2009) in CO J = 1 — 0, 
who used a similar SFR tracer. 

4.6. Systematics Affecting the Local SF Law 
4.6.1. Effect of the N on- detections 

In this study we analyze regions that have values over 
the adopted thresholds in both the Ssfr and £h 2 maps. 
To check for the effect of not including pixels that are 
detected in one axis but not the other, we include every 
pixel out to R25. This results in about 20 additional 
points, all with measurable SFR but no CO detection 
and thus having only upper limits for their gas surface 
densities. We find that these points closely follow the 
original distribution of points detected in both Esfr and 
Eh 2 at the limiting end of gas surface density. Thus, 
there are no new data trends hidden in the limits. They 
comprise only 3% of the total number of points obtained 
with the data selection criteria as mentioned in appendix 
[Til The impact of these points on the determination of 
the functional form of the SF law is negligible. 

4.6.2. Variations in the Data Selection Cuts 



It is a common practice to adopt one specific sensitivity 
limit in analyzing the gas-SFR surface density relation 
(for example, 3cr in Kennicutt et al. 2007; 2.5<r in Bigiel 
et al. 2008; 2a in Verley et al. 2010). The reason for 
adopting these sensitivity cuts is to ensure the reliability 
of the data. The choice of sensitivity limit, however, 
may have an impact on the determination of the local 
SF law, particularly given the limited dynamic range of 
the data. To explore the effect of this choice we have 
analyzed the gas-SFR surface density relation at several 
thresholds above 2a sensitivity. We find that the choice 
of limit has a measurable effect on the slope of the SF law 
for the pixel analysis, such that lower thresholds steepen 
the slope by as much as 30% — 40% (depending on the 
SFR tracer considered), with a simultaneous increase in 
the dispersion of the low surface brightness points. The 
determination of the slope for the aperture sampling and 
the azimuthally averaged radial profile are, on the other 
hand, robust to the choice of the sensitivity limit. 

4.6.3. Sensitivity to the Error Maps 

The results presented in this section are computed for 
a set of measurement error maps of SFR and gas surface 
densities. These maps are constructed under certain as- 
sumptions of the observational uncertainties. However, 
measurement uncertainties in flux calibration, continuum 
subtraction, and other parameters are propagated into 
the error maps. Variations in the assumptions made to 
include their contributions lead to changes in the error 
maps, which directly influence the regression analysis. 
For several sets of error maps with varying assumptions 
about the measurement uncertainties we find up to 40% 
variations in the slope measurements provided by the bi- 
variate regression methods. 

5. COMPARISON WITH PREVIOUS STUDIES 

In this section, we compare our results with recent 
studies of the spatially resolved SF law in nearby galaxies 
by Kennicutt et al. (2007), Bigiel et al. (2008), Blanc et 
al. (2009), and Verley et al. (2010). While making com- 
parison it should be borne in mind that, for a given SFR 
tracer and at a given kernel size, the table provides fitting 
results from three different bivariate regression methods. 
Our main results presented in various panels in Fig. 9, 
on the other hand, show the mean (Nmoi^t) and the dis- 
persion (<7 mo i,fit) of these three measurements. 
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Kennicutt et al. (2007) obtain a super-linear power law 
(N mo i ~ 1.37 ± 0.03) and an observed scatter of a- x ~ 0.4 
dex for NGC 5194 (M 51) using Ha + 24 fim as the SFR 
tracer, and apertures 520 pc in diameter centered on arm 
and inter- arm star-forming regions of M 51. They find a 
somewhat shallower power-law index in larger apertures 
(1850 pc in diameter). The authors subtract the diffuse 
component contribution using measurements in rectan- 
gular regions of the image, and employ FITEXY for the 
fitting. Although our aperture placement method is not 
strictly identical since we place apertures on the high 
surface density regions of NGC 4254, we find an approx- 
imately linear power law (Nmo^fit ± fmoi.fit ~ 0.9 ± 0.1) 
for our apertures independently of the DE subtraction in 
the SFR or molecular gas axis (see panel f in Fig. [9]). 
We also find that the functional form of the SF law is 
insensitive to the change of the aperture size, at least 
for the selected apertures in NGC 4254. The power law 
indices found at 500 pc and 1 kpc diameter apertures 
are well within their respective error bars (see panel f of 
Fig. [9]). The observed scatter is systematically smaller 
{o"i fit ~ 0.1 — 0.3 dex) in our study (see panel f of Fig. 
[HI. 

Bigiel et al. (2008) use the FUV + 24 ^m tracer, sim- 
ilar methodology in data analysis, OLS bisector fitting, 
and the same conversion factor, Xco- The authors find 
an approximately linear form (N mo i ~ 0.96±0.07) for the 
resolved molecular gas SF law in a sample of star-forming 
disk galaxies at 700 pc resolution. They measure a typi- 
cal molecular gas depletion timescale ~ 2 Gyr, under the 
assumption that the DIG contribution is negligible in 
both the FUV and MIR maps. They find a typical scat- 
ter of ~ 0.2 dex in the gas-SFR density relation. The 
authors did not consider the role of diffuse fraction in 
the SF law, therefore, their analysis is similar to the case 
/de = studied in this paper. Comparing our results of 
pixel analysis at /de = 0, we find a very similar power 
law index (N moUt ± er moljfit - 0.91 ± 0.03) in NGC 4254 
(see panel a of Fig. [9]). The molecular gas depletion 
timescale for this galaxy is T^ep ~ 1-7 Gyr. The estimate 
of intrinsic scatter obtained from OLS bisector method 
(<7i ~ 0.2 dex) is remarkably consistent with the Bigiel 
et al. (2008) study. 

Blanc et al. (2009) studied the central - 4 kpc of M 51 
using spectroscopic observations to obtain Shq- They 
find a slightly sub-linear (N mo i ~ 0.82 ± 0.05) functional 
form of the SF law with an intrinsic scatter u\ ~ 0.43 
dex. They place apertures 170 pc in diameter and mea- 
sure the DIG contribution using the [S II] /Ha ratio, 
finding /de ~ 11%. Despite the different DE removal 
and sampling _strategies, we also find an approximately 
linear slope (N mo i fit ± c mo i fit ~ 0.95 ± 0.2) at a simi- 
lar DIG fraction in NGC 4254 (see panel d of Fig. 0). 
The authors use a Monte Carlo method for fitting the 
SF law, which treats the intrinsic scatter in the relation 
as a free parameter. This method allows the inclusion of 
non-detections in both the Eh 2 and Ssfr maps, fitting 
the data in linear space. The advantage of this approach 
is that it is free from the systematics involved in per- 
forming linear regressions over incomplete data sets in 
logarithmic space. A drawback of this method is that it 
does not treat the data symmetrically, as it employs £h 2 
as the independent variable in the fits. Bivariate statis- 
tical methods such as the ones in our study do not easily 



allow for the inclusion of upper limits. Nonetheless, they 
permit a robust parametrization of the data with a sym- 
metric treatment of both axes, assuming that there is a 
good understanding of the data uncertainties. 

Verley et al. (2010) study the local SF law in M 33 
using azimuthally-averagcd radial profiles at a resolution 
of 240 pc, and using apertures sampling at various spa- 
tial scales (180 — 1440 pc). These authors also did not 
consider the role of diffuse fraction in the SF law. They 
use two different fitting techniques, including FITEXY, 
and several SFR tracers, including extinction-corrected 
Ha, extinction-corrected FUV, and a combination of ob- 
served FUV and total infrared luminosities. From aper- 
ture analysis they find that the molecular gas SF law is 
always super-linear and steeper (N mo i ~ 1.4—1.8) in Ha 
than in the other SFR tracers (N mol ~ 1.2 - 1.4). They 
find the same trend at all spatial scales. Their radial 
profile analysis, on the other hand, shows a shallower 
power-law index (N mo i ~ 1.1 — 1.3) by all three tracers. 
The results of our study are qualitatively similar to those 
of Verley et al. (2010). For example, our study also sug- 
gests that the power-law index is systematically steeper 
in Ha. We also obtain a shallower SF law when using 
azimuthally-averaged radial profiles. At coarser spatial 
resolutions, Verley et al. find a slight steepening of the 
SF law in all SFR tracers. However, we find little change 
in the power-law index over the explored range of spatial 
scales. 

6. SUMMARY AND CONCLUSIONS 

We study the spatially resolved molecular gas SF law 
in NGC 4254 within the central ~ 6.5 kpc, in an at- 
tempt to understand the combined effects of the under- 
lying assumptions on the functional form of the relation 
between molecular gas and SF. To this end, we use four 
different SFR tracers. These are Ha with azimuthally- 
averaged extinction correction, 24 fim, combined Ha and 
24 /zm, and combined FUV and 24 /im. We utilize vari- 
ous fitting procedures (the OLS bisector, FITEXY, and 
LINEXERR; described in appendix |E|) to constrain the 
parameters of the local SF law. We explore the effects 
of error weighting and signal-to-noise cuts on the results. 
We employ three different sampling strategies (pixel-by- 
pixel, aperture, and azimuthal averages) to probe the 
gas-SFR surface density relation at various spatial reso- 
lutions (500 pc and 1 kpc) and surface densities. 

We investigate the effect of diffuse emission on our abil- 
ity to measure the local SF law. Diffuse emission is an 
ubiquitous component of the maps in FUV, Ha, and 24 
/im, which may not be associated with star formation 
and comprises an unknown but perhaps significant frac- 
tion of the total disk emission (in the Milky Way, the 
diffuse emission in Ha related to the DIG constitutes 
10 — 15% of the total emission; Reynolds 1991). Simi- 
larly, there may be an analogous component of DE in the 
molecular gas axis of the SF law, comprised of unresolved 
clouds that are too small to host massive star formation. 
The contribution from DE is most important in the low 
surface brightness regime. To extract the DE from the 
maps we use spatial filtering (unsharp masking). 

We study the gas-SFR surface density relation for the 
molecular gas surface density range ~ 10 — 245 M pc~ 2 . 
This range is typical of normal star-forming galaxies but 
significantly smaller than starburst galaxies. The lower 
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limit of molecular gas surface density is consistent with 
recent studies which suggest that the atomic to molecular 
phase-transition occurs in the ISM at surface densities 
£ gas ~ 10 Mq pc -2 (Wong & Blitz 2002; Kennicutt et 
al. 2007; Leroy et al. 2008) 

Without additional data (for example, optical spec- 
troscopy) we cannot establish the fraction of DE cor- 
responding to the different SFR tracers in NGC 4254. 
Therefore, we take /de as an independent parameter and 
we explore the molecular gas-SFR surface density rela- 
tion for varying diffuse fraction. We find that in the high 
surface brightness regions sampled by our aperture anal- 
ysis (and dominating our azimuthal averages) the value 
of /de has little or no impact on the power-law index of 
the SF law, which is approximately linear in NGC 4254 
for all the SFR tracers considered (N^fit ~ 0.8 — 1.1; 
Fig. |9] right panels). Thus, the value of the logarithmic 
slope of the SF law can be robustly established in the 
high surface brightness regions of this galaxy. 

When lower surface brightness regions are included in 
the analysis, the removal of the DE, the choice of SFR 
tracer, and the choice of regression analysis all have an 
impact on the determination of the SF law (Fig. |9] left 
panels). For sub-dominant DE fractions in the SFR maps 
(/de 0-3), almost all tracers agree on an approximately 
linear SF law (Nmoi^t ~ 0.8—1.1). Removal of a DE com- 
ponent in both the SFR and the CO maps pushes this 
linear regime to slightly larger DE fractions (/de ^ 0.4; 
Fig. fl2|) . When the DE is assumed to dominate the 
emission (/de ^ 0.5) all tracers yield a super-linear SF 
law (Nmo^fit ^ 1.4). The discrepant tracer at low /de 
is Ha with a radially-dependent extinction correction. 
Since this is the SFR tracer that yields that largest scat- 
ter (<Xi,fit ~ 0.4 — 0.6 dex) in the gas-SFR surface den- 
sity relation, we think this steepness is an artifact of the 
combined effects of noise and performing the regression 
in logarithmic space. 

We explore a small range in spatial scales, as set by 
our choice of aperture sizes (500 pc and 1 kpc). In this 
range, we find that the results we obtain are independent 
of the spatial scale, with a very slight tendency to recover 
a steeper (but still approximately linear) SF law with a 
smaller scatter for larger aperture sizes (Fig. [9] right 
panels, points and dash-dotted lines). On the scales of 
the azimuthally-averaged radial profiles we recover an 
approximately linear SF law for all /de (N mo i fit ~ 1 — 1-2; 
Fig. [TO]). 

In NGC 4254 the intrinsic scatter (a - ;) of the SF law 
varies with the choice of the SFR tracer (Fig. ITT]) . In par- 
ticular, the 24 fim emission shows the tightest correlation 
with the molecular gas surface density (<Xi,fit ~ 0.1 — 0.4 
dex). This is likely because 24 /j,m emission is closely 
correlated with embedded SF, still associated with the 
parent molecular material. This suggests that 24 /jm is 
a very good tracer of SFR on timescales similar to the 
lifetimes of GMCs over spatial scales of several hundred 
parsecs. The combined Ha + 24 /jm tracer, which essen- 
tially applies a local extinction correction to Ha, yields 
the second tightest correlation (cxi,fit ~ 0.1 — 0.5 dex). 

The molecular SFE in this galaxy is constant out to 
R r-j O.5R25 (Fig. [13]). The disk averaged depletion 
timescale of the molecular gas (~ 2 Gyr) observed in 
NGC 4254 is fairly typical of large spirals (Bigiel et al. 



2008; Leroy et al. 2008) and in good agreement with 
previous observations (Wilson et al. 2009). The differ- 
ent SFR tracers agree very well with each other except 
for Ha, which yields a somewhat longer molecular gas 
depletion time. Like the power-law index, the SFE is 
independent of /de when the DE is sub-dominant. Re- 
moving a diffuse component from both the SFR and the 
molecular gas yields a SFE that is independent of /de- 
Although the presence of diffuse emission not associ- 
ated with star formation in the tracer used to determine 
the SFR or the molecular surface density should be a 
concern, at least in NGC 4254 the SF law can be de- 
termined in a precise and robust manner with no exact 
knowledge of the diffuse fraction in two cases: 1) in the 
high surface brightness regions independent of the value 
of /de, and 2) throughout the disk if /de < 30%. In 
both cases the resolved SF law is approximately linear in 
this galaxy. 
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Figure 2. Multi-wavelength maps of NGC 4254. The panels show FUV, Ha, and MIR 24 /im emission maps from left to right. The 
images are shown in logarithmic scale. The black contours represent the layout of CO J = 1 — map used in this study. The contours are 
in linear scale with the levels at 3 and 12 Jy Beam -1 km sec -1 . The figure shows a striking similarity between the distributions of hot 
dust traced by the MIR map and the cold molecular gas traced by CO J = 1 — in both upper an lower arm of the galaxy. 
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Figure 3. Two-dimensional unsharp masking. The top panels show the original FUV, Ha and MIR 24 (im emission maps from left to 
right. The middle panels show the corresponding smoothed images where the smoothing is performed by a 2D box kernel with the kernel 
size of 105" . The smoothed images are subtracted from the corresponding un-smoothed original images to produce unsharp masked maps 
(bottom panels). 
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Figure 4. Azimuthally averaged radial profiles of surface density maps shown at la up and below the mean trend at each radial point. 
Left Panel: The radial distribution of molecular gas surface density is shown for the most sensitive CARMA interferometer CO J = 1 — 
observation. The profiles derived from interferometer data alone as well as combining interferometer and single dish data are shown by 
gray shade and thin hatch, respectively. Right Panel: The radial profiles derived from Spuv+24/im, Sh , SH Q +24^m, and T,24fim maps 
are shown by long dashed line, thick hatched region, thin hatched region, and short dashed line, respectively. 
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Figure 6. Pixel-by-pixel analysis all four SFR tracers. The figure shows the molecular gas-SFR surface density relations at two represen- 
tative filter scales, 105" (filled contours) and 180" (black unfilled contours) as well as the case of no filtering, i.e., when /de=0 (colored 
unfilled contours). At any smoothing scale the amount of DE varies in the maps as shown at the bottom right in each panel (see also Table 
2). The gray scale represent the two dimensional histogram of the frequency of points, and the contours are placed at 90%, 75%, 50%, 
and 25% of the maximum frequency. The diagonal dotted lines represent constant SFE (e) where e ~100%, 10%, and 1% correspond to 
gas depletion time scale of 0.1, 1.0, and 10 Gyr. The filled circle represents the disk averaged estimate within the R25. The correlation 
strength varies among the maps. For a given SFR map, it depends on the degree of unsharp masking. The Pearson correlation coefficients 
at 105" and 180" are r ~ 0.4 — 0.55 (weak correlation) and ~ 0.55 — 0.7 (strong correlation), respectively, indicating how the subtraction of 
DE influence the correlation between the points. A comparison of the distributions of points shown by filled and unfilled contours clearly 
shows that the subtraction of DE affects the low surface density regions. 
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Figure 7. Aperture analysis for all SFR tracers for the selected apertures as shown in Fig. [5] The faint lines correspond to the contours 
of pixel distributions at 105" filter scale (shown as filled contours in Fig. [6}. The diagonal dotted lines present constant efficiencies as in 
Fig. [6] 
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Figure 8. Measurements provided by bivariate regression methods for SHn+24(im - Sh 2 relation using pixel-by-pixel sampling. The 
SFR tracer is subject to unsharp making with a smoothing scale of 105" which results in /dig = 0-45 and /mir=0.56. The gray scale 
represent the two dimensional histogram of the frequency of points, and the contours are placed at 90%, 75%, 50%, and 25% of the 
maximum frequency. For this particular distribution of points, the FITEXY provides the shallowest slope (N mo i ~ 1.3 short dash line) 
The OLS bisector (N mo i ~ 1.6 dashed-dot line) and the LINEXERR (N mo i ~ 1.4 long dash line) methods, on the other hand, provide 
steeper estimates. The estimated slope, zero point, and formal error provided by individual fit arc shown at the lower right corner. The 
intrinsic scatter provided by FITEXY, OLS bisector, and LINEXERR methods are a\ ~ 0.4 0.3, and 0.3 dex, respectively. The vertical 
and horizontal axes range from ~ 0.1 — 114 Mq Gyr -1 pc~ 2 and ~ 10 — 245 Mq pc -2 , respectively. 
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Figure 9. Dependence of the power law index (N mo i) of the molecular gas SF law on various diffuse fractions. The figure summarizes 
results for both pixel (left panels) and aperture (right panels) samplings. Results from 6" ~ 500 pc (points) and 12" ~ 1 kpc (dashed-dot 
line) diameter aperture samplings are shown in the right panels. The methodological mean (N mo i fj t ) and dispersion (<r mol fj t ) from three 
measurements are shown, respectively, by gray filled circles and vertical lines. The black points connected by dashed lines in panels a 
and b (e and f) represent the locus of N mol fj t at /fuv(/dig)- The horizontal dotted lines represent N mo i = 1 and N mol = 1.5. The 
vertical dotted lines demarcate the regions where the DE is assumed to be (from left to right) the dominant, significant, sub-dominant, and 
negligible component of the total disk emission. Each point along the horizontal axis has 15% uncertainty (see appendix [F] for details). 
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Figure 10. Azimuthally averaged radial profile analysis all four SFR tracers. Figure shows the radial relations at the dominant (filled 
circle), significant (open square) and negligible (/de = 0) diffuse fractions. The faint lines correspond to the contours of pixel distributions 
at 105" filter scale (shown as filled contours in Fig. |BJ. The diagonal dotted lines represent constant efficiencies as in Fig. [6] The bin size 
is 500 pc. The absence of an extended tail along the vertical axis at the low surface density regions, irrespective of the subtraction of the 
amount of diffuse fractions, indicates that azimuthally averaged radial profile systematically sample the high SFR surface density regions 
compared to pixel sampling. See Table [7] for the derived fitted parameters for all SFR tracer maps at these three diffuse fractions. 
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Figure 11. The intrinsic scatter (<Ti) in the gas-SFR surface density relations at various diffuse fractions for pixel (left panels) and 
aperture (right panels) samplings. Results from 6" ~ 500 pc (points) and 12" ~ 1 kpc (dashed-dot line) diameter aperture samplings are 
shown in the right panels. Figure shows the methodological mean of the observed scatter (<7j a t ; "mean dispersion") by filled circle and the 
dispersion in the observed scatter ("dispersion of dispersion") from linear regression methods by vertical solid lines. Figure highlights that 
the estimate of intrinsic scatter depends on the choice of the SFR tracer as well as on the treatment of DE. For example, the scatter is 
systematically larger in the Sjj 2 - relation (methodological dispersion, ctj fi t ~ 0.3 — 0.6 dex). The scatter becomes smaller for larger 

aperture (dashed-dot line in the right panel) because of the reduction in the sampled points. The vertical lines have similar meanings as 
in Fig. \9\ Each point along the horizontal axis has 15% uncertainty (see appendix [F] for details). 
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Figure 12. Dependence of the power law index (N mo i) of molecular gas SF law when both SFR tracer and CO J = 1 — map are subject 
to unsharp masking. The dashed and thin solid lines, respectively, show the methodological dispersions when the selection threshold for 
the data in the Ejj 2 unsharp-maskcd map is placed at lcr (3.7 Mq pc -2 ) and 2.5<r (9.2 Mq pc — 2 ). The thick solid line illustrates the case 
when only the SFR tracers are subject to masking. Presentation style is similar to Fig. [9] 
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Figure 13. Disk averaged SFE (Gyr 1 ) and the depletion timescale, r^ep (Gyr), at various diffuse fractions. The SFE and the r^ &v are 
shown on the left and right vertical axis, respectively. Results for both (6" by points and 12" by dashed-dot line) aperture samplings are 
shown in the right panels. Parameters estimated from applying unsharp masking to the SFR tracer only, and to both the SFR and the gas 
maps, are shown in black and gray respectively. For a given SFR tracer, the subtraction of DE introduces ~50% (Sh q ) up to a factor 3 
variations in the global SFE. At a given /de, the global SFE derived from four tracers are approximately consistent with one another (see 
Table E), although Ha yields a lower depletion time. The dashed lines added by black points shown in Figs. f9l and fill have omitted for 
clarity of presentation. The vertical lines have similar meanings as in Fig. [9] 
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Table 4 

Fitted Parameters from Pixel Analysis 



Filt. I /dig FUV+24 )im Ha Ha+24 )im 24 |im 



(") 




Log A 


N mo[ 




Log A 


N mo , 


(7 \ 


Log A 


N mo i 




Log A 


N mo[ 


(Ti 


75 


0.55 


-1.97 ± 0.11 
—0.96 ± 0.14 
-1.36 ± 0.16 


1.71 ± 0.06 
1.15 ± 0.07 
1.35 ± 0.09 


0.50 
0.37 
0.36 


-1.99 ± 0.13 
— 0.83 ± 0.20 
-1.29 ± 0.21 


1.67 ± 0.07 
1.02 ± 0.11 
1.25 ± 0.11 


0.54 
0.47 
0.47 


-2.04 ± 0.11 
— 1.12 ± 0.17 
-1.55 ± 0.18 


1.70 ± 0.06 
1.19 ± 0.09 
1.40 ± 0.10 


0.49 
0.40 
0.41 


-1.75 ± 0.10 
— 1.00 ± 0.13 
-1.40 ± 0.15 


1.57 ± 0.05 
1.15 ± 0.07 
1.35 ± 0.08 


0.42 
0.33 
0.31 


105 


0.45 


-1.57 ± 0.10 
-1.06 ± 0.13 
-1.26 ± 0.13 


1.57 ± 0.05 
1.27 ± 0.07 
1.36 ± 0.07 


0.40 
0.30 
0.30 


-1.79 ± 0.10 
-0.77 ± 0.16 
-1.22 ± 0.18 


1.60 ± 0.05 
1.02 ± 0.09 
1.25 ± 0.10 


0.50 
0.42 
0.42 


-1.66 ± 0.09 
-1.13 ± 0.15 
-1.45 ± 0.14 


1.57 ± 0.05 
1.26 ± 0.08 
1.41 ± 0.08 


0.41 
0.32 
0.33 


-1.27 ± 0.08 
-0.95 ± 0.12 
-1.19 ± 0.11 


1.38 ± 0.04 
1.19 ± 0.06 
1.30 ± 0.06 


0.33 
0.23 
0.23 


180 


0.21 


-0.50 ± 0.05 
-0.44 ± 0.09 
-0.58 ± 0.08 


1.10 ± 0.03 
1.05 ± 0.05 

1.11 ± 0.04 


0.24 
0.15 
0.15 


-1.28 ± 0.07 
-0.83 ± 0.13 
-1.10 ± 0.12 


1.41 ± 0.04 
1.15 ± 0.07 
1.27 ± 0.07 


0.36 
0.30 
0.30 


-0.68 ± 0.05 
-0.57 ± 0.09 
-0.74 ± 0.08 


1.15 ± 0.05 
1.07 ± 0.05 
1.15 ± 0.04 


0.25 
0.17 
0.17 


-0.34 ± 0.04 
-0.42 ± 0.07 
-0.53 ± 0.07 


0.98 ± 0.02 
1.00 ±0.04 
1.05 ± 0.04 


0.20 
0.10 
0.10 




0.0 


-0.03 ± 0.04 
-0.06 ±0.06 
-0.16 ± 0.06 


0.90 ± 0.02 
0.90 ± 0.03 
0.94 ±0.03 


0.18 
0.10 
0.10 


-1.10 ± 0.06 
-0.74 ± 0.11 
-0.97± 0.11 


1.33 ± 0.03 
1.12 ± 0.06 
1.23 ± 0.06 


0.33 
0.27 
0.26 


-0.36 ± 0.04 
-0.34 ± 0.07 
-0.47 ± 0.07 


1.01 ± 0.02 
0.98 ± 0.04 
1.04 ± 0.04 


0.21 
0.13 
0.13 


-0.02 ± 0.04 
-0.07 ±0.05 
-0.16 ± 0.05 


0.82 ± 0.02 
0.86 ± 0.03 
0.90 ± 0.03 


0.16 
0.10 
0.10 



Note. — The normalization constant, power-law index, and intrinsic scatter of the molecular gas SF law are derived from unsharp masked 
SFR tracers (i.e., unsharp masking along the vertical axis). The filter scales from the top correspond to the dominant (/dig ^ 50%), significant 
(/dig ~ 30 — 50%), and sub-dominant (/dig ~ 10 — 30%) diffuse fractions. The bottom row shown by the dashed line represents /dig— 0. At 
a given filter scale, the three rows under FUV+24/i.m provide the measurements of the OLS bisector, FITEXY, and LINEXERR fitting method, 
respectively. Typical uncertainty in the intrinsic scatter is few percent. 



Table 5 

Fitted Parameters from 6" Aperture Analysis 



Filt. I /dig FUV+24 ^m Ha Ha+24 ^m 24 ^m 

(") Log A N mol trj Log A N mo i erj Log A N mol trj Log A N mol <rT 



75 


0.55 


-0.09 ±0.17 1.29 ± 0.09 0.28 
-0.26 ± 0.23 0.84 ±0.12 0.25 
-0.30 ± 0.26 0.86 ±0.14 0.23 


-1.61 ± 0.22 1.49 ±0.11 0.37 
-0.40 ± 0.37 0.85 ± 0.20 033 
-0.51 ± 0.35 0.90 ±0.19 0.33 


-1.27 ±0.17 1.34 ± 0.08 0.30 
-0.30 ± 0.23 0.82 ±0.12 0.20 
-0.39 ± 0.28 0.86 ±0.15 0.30 


-1.03 + 0.16 1.23 + 0.08 0.25 
-0.24 + 0.23 0.80 + 0.12 0.20 
-0.35 + 0.24 0.86 + 0.13 0.20 


105 


0.45 


-0.67 ±0.14 1.13 ± 0.07 0.22 
-0.03 ± 0.16 0.78 ±0.08 0.16 
-0.13 ± 0.20 0.84 ±0.11 0.17 


-0.34 ±0.18 1.40 ± 0.09 0.32 
-0.33 ± 0.30 0.86 ±0.16 0.27 
-0.42 ± 0.30 0.90 ±0.16 0.28 


-0.84 ±0.14 1.17 ± 0.07 0.23 
-0.10 ±0.21 0.78 ±0.11 0.18 
-0.19 ± 0.22 0.82 ±0.12 0.19 


-0.61 + 0.13 1.06 + 0.07 0.20 
-0.06 + 0.18 0.77 + 0.10 0.16 
-0.14 + 0.19 0.81 + 0.10 0.14 


180 


0.21 


-0.25 ± 0.12 0.99 ±0.06 0.17 
+0.12 ± 0.16 0.79 ±0.08 0.12 
-0.04 ± 0.16 0.82 ±0.08 0.11 


-0.87 ±0.14 1.22 ± 0.07 0.24 
-0.16 ± 0.20 0.84 ±0.11 0.22 
-0.23 ± 0.23 0.87 ±0.12 0.20 


-0.35 ±0.11 0.99 ± 0.06 0.17 
+0.07 + 0.15 0.76 + 0.08 0.12 
-0.02 + 0.17 0.79 + 0.09 0.12 


-0.20 + 0.10 0.92 + 0.05 0.15 
+0.06 + 0.15 0.78 + 0.08 0.10 
-0.05 + 0.14 0.78 + 0.07 0.10 




0.0 


-0.03 ± 0.10 0.88 ± 0.05 0.14 
+0.29 ±0.12 0.74 ± 0.07 0.10 
+0.24 + 0.14 0.76 + 0.07 0.10 


-0.83 ±0.15 1.21 ± 0.08 0.23 
-0.20 ± 0.18 0.88 ± 0.09 0.20 
-0.25 ± 0.22 0.90 ±0.12 0.20 


-0.21 + 0.11 0.95 + 0.06 0.16 
+0.16 + 0.16 0.75 + 0.08 0.10 
-0.08 + 0.15 0.78 + 0.08 0.10 


+0.03 + 0.09 0.83 + 0.05 0.13 
+0.24 + 0.13 0.71 + 0.06 0.10 
+0.22 + 0.12 0.72 + 0.06 0.10 



Note. — The normalization constant, power-law index, and intrinsic scatter of the molecular gas SF law are derived from unsharp masked 
SFR tracers (i.e., unsharp masking along the vertical axis). The filter scales from the top correspond to the dominant (/dig ^ 50%), significant 
(/dig ~ 30 — 50%), and sub-dominant (/dig ~ 10 — 30%) diffuse fractions. The bottom row shown by the dashed line represents /dig— 0. At 
a given filter scale, the three rows under FUV+24/j.m provide the measurements of the OLS bisector, FITEXY, and LINEXERR fitting method, 
respectively. Typical uncertainty in the intrinsic scatter is few percent. 
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Table 6 

Fitted Parameters from 12" Aperture Analysis 



Filt. I /dig FUV+24 )im Ha Ha+24 )im 24 fim 

(") Log A N mo i ct! Log A N moI Log A N mo i ctI Log A N mol 57 



75 


0.55 


-1.20 ± 0.30 1.33 ±0.15 0.23 
-0.12 ± 0.44 0.76 ± 0.23 0.20 
-0.10 ± 0.50 0.75 ± 0.26 0.20 


-1.42 ± 0.33 1.39 ±0.17 0.26 
-0.31 ± 0.67 0.81 ± 0.36 0.20 
-0.14 ± 0.56 0.73 ± 0.29 0.25 


-1.30 ± 0.33 1.34 ±0.16 0.26 
-0.05 ± 0.52 0.63 ± 0.28 0.20 
+0.04 + 0.56 0.64 ± 0.30 0.24 


-1.22 ±0.30 1.31 ±0.15 0.25 
-0.11 ± 0.50 0.73 ± 0.26 0.20 
-0.07 ± 0.49 0.71 ± 0.26 0.20 


105 


0.45 


-0.88 ± 0.27 1.22 ±0.13 0.18 
-0.39 ± 0.84 0.99 ±0.42 0.15 
-0.15 ± 0.39 0.84 ±0.21 0.15 


-1.45 ± 0.35 1.44 ±0.18 0.25 
-0.21 ± 0.85 0.82 ± 0.43 0.15 
-0.22 ± 0.56 0.80 ± 0.29 0.23 


-1.01 ± 0.26 1.25 ±0.13 0.20 
+0.25 + 0.73 0.62 ± 0.37 0.10 
-0.02 ± 0.44 0.72 ± 0.23 0.18 


-0.84 ±0.25 1.17 ±0.12 0.21 
+0.30 + 0.42 0.61 +0.21 0.10 
-0.09 + 0.37 0.77 + 0.20 0.14 


180 


0.21 


-0.41 ±0.21 1.06 ±0.11 0.12 
-0.03 ± 0.22 0.82 ±0.12 0.10 
-0.05 ± 0.29 0.86 ±0.15 0.10 


-1.01 ± 0.25 1.28 ±0.13 0.20 
-0.17 ± 0.44 0.83 ± 0.23 0.12 
-0.16 ± 0.42 0.83 ± 0.22 0.17 


-0.54 ±0.21 1.08 ±0.11 0.15 
-0.09 ± 0.27 0.75 ±0.14 0.10 
-0.06 ± 0.32 0.76 ±0.17 0.11 


-0.39 + 0.21 1.01 + 0.10 0.12 
-0.07 + 0.25 0.84 + 0.14 0.10 
+0.02 + 0.27 0.79 + 0.14 0.10 




0.0 


-0.09 ± 0.18 0.93 ±0.09 0.10 
-0.31 ±0.31 0.71 ±0.17 0.10 
+0.17 ± 0.25 0.79 ±0.13 0.10 


-0.94 ± 0.26 1.26 ±0.13 0.18 
-0.16 ± 0.32 0.85 ±0.17 0.13 
-0.20 ± 0.39 0.87 ± 0.20 0.15 


-0.34 ± 0.20 1.00 ±0.10 0.12 
-0.19 ± 0.30 0.72 ±0.16 0.10 
-0.12 ± 0.28 0.76 ±0.15 0.10 


-0.11 + 0.18 0.90 + 0.01 0.10 
-0.15 + 0.15 0.76 + 0.08 0.10 
+0.21 + 0.24 0.73 + 0.13 0.10 



Note. — The normalization constant, power-law index, and intrinsic scatter of the molecular gas SF law are derived from unsharp masked 
SFR tracers (i.e., unsharp masking along the vertical axis). The filter scales from the top correspond to the dominant (/dig ^ 50%), significant 
(/dig ~ 30 — 50%), and sub-dominant (/dig ~ 10 — 30%) diffuse fractions. The bottom row shown by the dashed line represents /dig— 0. At 
a given filter scale, the three rows under FUV-|-24/j,m provide the measurements of the OLS bisector, FITEXY, and LINEXERR fitting method, 
respectively. Typical uncertainty in the intrinsic scatter is few percent. 



Table 7 

Fitted Parameters from Azimuthally Averaged Radial Profile Analysis 



Filt. 


/dig 


FUV+24 fj,m 




Ha 




Ha+24 (im 




24 fj,m 




(") 




Log A 


N mol 


Log A 


N mo i 


Log A 


N mol 


Log A 


N mo i 


75 
105 
180 


0.55 
0.45 
0.21 
0.0 


-0.82 ± 0.11 
-0.78 ± 0.08 
-0.53 ± 0.04 
-0.16 ± 0.03 


1.13 ± 0.07 
1.17 ± 0.05 
1.12 ± 0.02 
0.98 ± 0.02 


-0.81 ± 0.15 
-0.78 ± 0.12 
-0.63 ± 0.08 
-0.64 ± 0.06 


1.09 ± 0.09 

1.11 ± 0.07 
1.09 ± 0.05 

1.12 ± 0.04 


-0.79 ± 0.13 
-0.76 ± 0.10 
-0.55 ± 0.05 
-0.37 ± 0.04 


1.07 ± 0.08 
1.11 ± 0.06 

1.08 ± 0.03 
1.02 ± 0.02 


-0.94 ± 0.11 
-0.87 ± 0.07 
-0.57 ± 0.03 
-0.20 ± 0.02 


1.16 ± 0.06 
1.18 ± 0.04 
1.10 ± 0.02 
0.95 ± 0.01 


75 
105 
180 


0.55 
0.45 
0.21 
0.0 


-0.16 ± 0.03 
-0.35 ± 0.06 
-0.20 ± 0.10 
-0.44 ± 0.08 


0.98 ± 0.02 

1.07 ± 0.04 

1.08 ± 0.06 
1.12 ± 0.04 


-0.64 ± 0.06 
-0.38 ± 0.08 
-0.33 ± 0.08 
-0.48 ± 0.08 


1.12 ± 0.04 
1.04 ±0.05 

1.07 ± 0.05 

1.08 ± 0.05 


-0.37 ± 0.04 
-0.35 ± 0.06 
-0.24 ± 0.08 
-0.46 ± 0.06 


1.02 ± 0.02 
1.02 ± 0.04 
1.05 ± 0.06 
1.08 ± 0.04 


-0.20 ± 0.02 
-0.45 ± 0.05 
-0.24 ± 0.10 
-0.53 ± 0.07 


0.95 ± 0.01 
1.09 ± 0.03 
1.06 ± 0.04 
1.14 ± 0.04 



Note. — Top: parameters in the top four rows are derived from unsharp masked SFR tracer maps. Bottom: parameters in the bottom four rows 
are derived when both molecular gas and SFR maps arc subject to unsharp masking. The intrinsic scatter a- x is in radial profile analysis and hence 
the parameters are obtained from the OLS bisector method. The filter scales from the top correspond to the dominant (/dig ^ 50%), significant 
(/dig ~ 30 — 50%), and sub-dominant ((/dig ~ 10 — 30%) diffuse fractions. The bottom row shown by the dashed line represents /dig— 0. 
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APPENDIX 
A. NGC 4254 

The optical structure of NGC 4254 shows a one-armed appearance (m = 1 mode), unlike most grand-design spirals 
where symmetric modes are prominent. This unusual morphology has made it the target of several observational and 
numerical studies (Phookun et al. 1993; Chcmin et al. 2006; Sofue et al. 2003). Various explanations for its asymmetry 
have been put forth including the superposition of spiral modes induced by global gravitational instability (lye et 
al. 1982), the asymmetric accretion of gas onto the disk (Phookun et al. 1993; Bournaud et al. 2005), ram pressure 
stripping (Sofue et al. 2003; Kantharia et al. 2008), a close high-speed encounter plus ram pressure (Vollmer et al. 
2005), a high-speed encounter only (Due & Bournaud 2008), and harassment while entering Virgo (Haynes et al. 2007). 
On the other hand, detailed photometric studies of NGC 4254 show that it is photometrically similar to other Sc type 
spirals, and it has no close companion (Phookun et al. 1993). 

NGC 4254 has a flat rotation curve with v rot ~ 150 km sec -1 up to 200" (~ 16 kpc) from its center (Guhathakurta 
et al. 1988). Apart from its morphological peculiarity, it does not show any anomaly in its properties. For example, 
the distribution of SF through out the disk based on Ha emission has been classified as normal (Koopman & Kenney 
2004). Its molecular gas fraction in the disk is similar to field galaxies, suggesting that it is entering Virgo for the first 
time and external agents (such as ram pressure) have not yet been effective at striping its gas, despite its morphology 
(Nakasihsi et al. 2006). NGC 4254 is a metal rich galaxy with a median 12+log[0/H] ~ 9.0 across the disk (Vila-Costas 
& Edmunds 1992; Zaritsky et al. 1994). 

B. SFR SURFACE DENSITY MAPS 

We construct four different SFR maps from combining FUV and MIR maps (Leroy et al. 2008), optical Ha emission 
map (Kennicutt et al. 2007; Prescott et al. 2007), optical Ha emission and MIR maps (Calzetti et al. 2007), and 
MIR map (Calzetti et al. 2007). The details of the construction of these SFR tracers can be found in the references. 
These different SFR tracers probe different time scales and hence the SF history of any particular galaxy. For 
example, the Ha emission traces gas ionized by massive (M > 10 M© ) stars over a timescale of < 20 Myr. The 
FUV luminosity corresponds relatively older (< 100 Myr), less massive (M > 5 M Q ) stellar populations. The MIR 
24 /im emission mostly traces re-processed radiation of newborn (few Myr) OB stellar associations embedded inside 
the parent molecular clouds. Although star clusters break from their parent clouds in less than 1 Myr, they remain 
associated with it for a much longer time scale, t <~ 10 — 30 Myr. Thus this is the time scale associated with the MIR 
24 /im emission as a tracer. However, if there is significant heating of the small dust grains by late B and A stars, the 
time scales associated with 24 /im emission may be as long as ~ 100 Myr. For composite SFR tracer such as FUV + 
24 /im or Ha + 24 /im while the MIR emission traces the dust-obscured fraction of the SF, and the un-obscured SF 
can be traced either by FUV continuum or optical nebular emission. 

Each pixel in the FUV, Ha, and 24 /im emission maps has some uncertainties coming from observation and data 
processing. These errors propagate to the uncertainties in the surface density maps. For example, uncertainties in the 
flux calibration, the stellar continuum subtraction, and the background subtraction contribute to the total error budget 
of Ha. We assume 10% uncertainty for each of these components and construct Ha error map ctho by adding these 
terms quadratically. A flat 10% uncertainty in stellar continuum subtraction for the entire galaxy disk is probably 
unrealistic because stellar emission dominates in the galaxy center. Therefore, it will have higher contribution at the 
high surface brightness regions. We take this limit on ad hoc basis. The assumed fraction in the continuum subtraction 
may a lower limit since it varies along the galaxy disk contributing as much as 30% to In a flux uncertainty (Koopman 
et al. 2001). Finally, there is a considerable (calibration) uncertainty in Ha to SFR conversion (Calzetti et al. 2007). 
We take all these factors into consideration when constructing the SFR error map (au a ) from the Ha image. The 
FUV and mid-IR error maps (ofuv an d 0mir) are constructed in a similar manner but without the contribution from 
the stellar continuum. The error in the composite SFR maps (FUV + 24 /im or Ha + 24 /im) are constructed by 
combining the appropriate terms from the respective images. 

C. MOLECULAR GAS SURFACE DENSITY MAP 

A conversion factor (Xco) is frequently used to determine the distribution of molecular hydrogen (H2) from the 
CO J = 1 — images. In this study, we use Xco = 2.0 x 10 20 to be consistent with several other current studies 
(Wong & Blitz 2002; Komugi et al. 2005; Gardan et al. 2007; Bigiel et al. 2008; Leroy et al. 2008). The choice of the 
calibration factor linearly scales the estimated gas densities. It is worthwhile asking how accurate is the assumption of 
a single conversion factor. Although the CO-to-H2 appears to be approximately constant for resolved GMCs (Bolatto 
et al. 2008), on the hundreds of parsecs scales sampled in this study it may vary across the disk, especially in the 
low density low metallicity regions (e.g., Garcia-Burillo et al. 1993). We adopt a constant Xco f° r simplicity and 
consistency with most recent studies (Wong & Blitz 2002; Komugi et al. 2005; Gardan et al. 2007; Bigiel et al. 2008; 
Leroy et al. 2008). Furthermore, we expect the effects of a varying conversion factor to be most important in the 
outer portions of the disk, which are beyond the central region sampled by this study. The H 2 surface densities are 
multiplied by a factor of 1.40 to account for the mass contribution of helium and are expressed in units of M Q pc -2 . 

To obtain the integrated (interferometer + single) CO map, first, the CARMA cube was de-convolved using the 
implementation of the CLEAN algorithm in the MIRIAD task mossdi and the velocity channels (2.6 km sec -1 ) in the 
IRAM cube has been re-binned to have the same velocity resolution (10 km sec -1 ) of the CARMA cube. The MIRIAD 
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task immerge was then used to combine the CARMA and the IRAM cube in the image plane (e.g., Stanimirovic et 
al. 1999). The integrated intensity map was created by convolving each plane of the original data-cube with a 15" 
Gaussian to degrade its resolution, selecting regions in each velocity plane where the signal is larger than 2.5a, and 
using these regions as masks on the original cube to compute an integrated intensity. The corresponding error map 
was computed by considering the la rms value of 22 mJy Beam -1 km sec -1 for each plane and multiplying by the 
square root of the number of planes that make up each integrated intensity point. Although the RMS in each plane is 
approximately constant within the field-of-view of the mosaic, the noise in the integrated intensity map is a function 
of position because of the change in line width of the CO emission. 

D. DATA SAMPLING 

Existing studies of spatially resolved SF law in galaxies use three different data sampling methodologies. These are: 
azimuthally averaged radial profile (Kennicutt 1989; Martin & Kennicutt 2001; Wong & Blitz 2002; Boissier et al. 2003; 
Heyer et al. 2004; Schuster et al. 2007; Bigiel et al. 2008); aperture analysis encompassing the star forming regions and 
centering on Ha and MIR 24 /im emission peaks (Kennicutt et al. 2007; Blanc et al. 2009); and finally, pixel-by-pixel 
analysis (Bigiel et al. 2008; Leroy et al. 2008). To understand the methodological impact on the determination of the 
local SF law we have incorporated all three methods in our study since each of these methodologies has strengths and 
weaknesses. 

Pixel analysis probes the gas-SFR density relation at the smallest spatial resolution. Aperture averages, on the 
other hand, provide information relevant to ensemble averages of representative regions. Azimuthally averaged radial 
profiles, in general, provide estimates over much larger scales. This method of sampling produces overall radial trends 
suppressing local variations. Radial profiles, by construction, are less susceptible to angular resolution. Although the 
first sampling probes the smallest spatial scales in the map, it is heavily dominated by its low surface brightness regions. 
The latter two samplings, by construction, probe only the high surface brightness regions where our signal-to-noise is 
best and DE in both the SFR and the molecular gas tracers is likely to be a minor contaminant. 

In our pixel analysis, we consider regions of the galaxy where both Esfr and Eh 2 are higher than 2.5a. All images 
are resolution- matched to 6", and re-sampled at the Nyquist rate (3" x 3" pixels). For the aperture analysis we lay 
down a total of 102 circular apertures 6" (500 pc) in diameter, in regions of high surface density, mostly following the 
spiral arms (see Fig. [5J- These apertures are approximately independent, with only a small fraction (5%) of overlap. 
We also consider 34 circular apertures 12" (1 kpc) in diameter with a similar overlapping fraction. It should be stressed 
here that sampling the same region in all four SFR surface density maps simultaneously determines the total number 
of apertures. The aperture centers for both small and large radii were chosen independently. For any Esfr map all the 
apertures combined cover approximately « 45% of the total area covered by all the pixels above 2.5er limit. Because of 
the nature of the selection method the aperture samplings contain ~ 65% of the total emission when compared pixel 
sampling. 

E. STATISTICAL METHODOLOGY AND FITTING PROCEDURES 

The SF law, as expressed by Eq. 1, is a power law relationship between the SFR surface density and the gas surface 
density. The outcome of any regression analysis, with the object of finding A and N, depends on the treatment of 
the data and its measurement errors, and on the intrinsic scatter of the observables. The intrinsic scatter reflects the 
variations of local physical properties of star-forming regions (for example, evolutionary stages, stellar populations, 
metallicity, obscuration, etc). The measurement errors, on the other hand, depend entirely upon the observations 
and the subsequent data reduction. The factors that contribute to the measurement errors include flux calibration, 
continuum subtraction, and background subtraction. 

The linear regression methods can be divided into two broad classes depending on whether the intrinsic scatter or 
the measurement error dominate. When seeking a linear relationship between two variables in a data set where both 
of them have small measurement errors but substantial yet unknown intrinsic scatter, the ordinary least square (OLS) 
bisector method provides one the best solutions, (Isobe et al. 1990; Feigelson & Babu 1992). 

When both the dependent and the independent variables are subject to measurement errors and intrinsic scatter 
of comparable magnitude, the regression analysis becomes more complex. Several bivariate regression methods have 
been developed to deal with astronomical problems but each method has its own advantages and disadvantages (see 
Feigelson & Babu 1992 for a detailed account). The most widely used bivariate regression analysis is based on the 
least squares technique (FITEXY; Press et al.1992 and references therein). Akritas & Bershady (1996) extended 
the OLS bisector method incorporating measurement error and the intrinsic scatter. This estimator is known as 
the bivariate correlated error and scatter (BCES) estimator. Kelly (2007) developed a bivariate estimator based on 
Bayesian statistics (LINEXERR). 

The FITEXY and LINEXERR estimators differ in the underlying assumption of the nature of the true relationship 
between independent and dependent variables. Both the BCES and LINEXERR estimators assume that the data 
points are scattered around the true linear relationship. The FITEXY, on the other hand, assumes that there is no 
intrinsic scatter and the true points lie exactly on a straight line, providing a solution for data showing a perfect 
correlation. While the BCES and LINEXERR estimators incorporate the correlated measurement error, FITEXY 
does not account for it. Correlated error arises when the dependent and the independent variables both arc subjected 
to the same uncertainty. For example, surface densities are defined as the ratio of the total mass to the de-projected 
area of the disk. An uncertainty in the inclination measurement leads to a correlated measurement error in the surface 
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density of gas and SFR. The covariance term broadens the actual distribution of data points and thus provides a flatter 
relationship if unaccounted for (Akritas & Bershady 1996). 

In this study we have used the OLS bisector, FITEXY, and LINEXERR estimators. We note here that while the 
OLS bisector method will guard our analysis against possible flaw in constructing measurement error maps, the other 
two methods will be required for a complete analysis of data. Our experience shows that the BCES bisector method is 
highly sensitive to the measurement error resulting in unstable slope estimates compared to the three other methods 
mentioned above. We do not use this estimator in our analysis. 

We construct surface density measurement error maps as mentioned in section ij2j The covariant term for each 
pixel is calculated from the error map. The intrinsic scatter is directly provided by the LINEXERR estimator. For 
the OLS bisector and FITEXY estimators we estimate it using the best fit line. To include scatter in the FITEXY 
estimator we iteratively adjust the error along the Y-axis iteratively until we achieve a reduced chi-squared ~ 1. For 
this estimator, the total dispersion along the y-axis is o~y — 4- of. The adjustment is made for o\ which is a measure 
of the intrinsic dispersion in the gas-SFR surface density relation, while a m is the measurement uncertainty obtained 
from error propagation. The formal error in each parameter quoted in this paper comes from the bootstrap sampling 
of 1000 realizations of the data points. 

F. UNSHARP MASKING 

We use unsharp masking to model and remove the local variations of DE from the maps. The choice of the size of 
the median filter kernel plays a vital role in selecting the amount of diffuse component in the total disk emission. As 
mentioned earlier, we explore a number of filter sizes in each SFR tracer, carrying out our analysis for each case (see 
Table [5]) . We should stress here that our intention is to use the simplest model for a reliable estimate of the diffuse 
fraction in the disk of NGC 4254 which can be used with a reasonable confidence to explore the main goal of this 
study. 

For a filter size of 75" the fraction of the total emission contained in the smooth Ha map is /dig ~ 60%, on the 
acceptable high-end of the DIG fractions observed in other galaxies (Kennicutt et al. 1995; Ferguson et al. 1996). 
Thus, DE in Ha is likely overestimated for kernel sizes below 75" (corresponding to physical scales ~ 6 kpc). Using 
smaller kernels for unsharp masking results in over-subtracted maps containing undesirable artifacts (many negative 
pixels and increased scatter). The upper end of acceptable kernel sizes is determined by the image size. We set it 
roughly to two-third the image size, corresponding to 210" or a physical length of ~ 13 kpc. With this scale we obtain 
/dig ~ 15% which is approximately the /dig observed in the Galaxy (Reynolds 1991). Similar kernel sizes are also 
explored in the FUV and 24 /jm maps. As expected, the diffuse fractions in the FUV, Ha, and 24 /im emission maps 
vary substantially from one another. Within our range of kernel sizes, /fuv ~ 40 — 75%, while /mir ~ 20 — 70%. 

We show our azimuthally averaged radial distributions of DE fractions for NGC 4254 in Fig. [T3J Panels A, B, and 
C show the radial profiles for /fuv, /dig, and /mir respectively. In each panel, we show our results for six different 
smoothing filters. 

It is evident from the figure that the /de profiles increase radially outward. In spite of their distinct origins, all three 
diffuse fractions show this remarkably similar trend. Beyond ~10 kpc these profiles flatten out. For smaller smoothing 
scales, /mir varies as much as a factor of four whereas /dig varies up to a factor three. The /fuv shows a factor of 
two radial variation. For longer smoothing scales the variation is about a factor of two for all /de- For /dig, a similar 
trend was observed in the Sculptor group spiral NGC 7793 (Ferguson et al. 1996). Panel D shows the average fraction 
of diffuse emission, Joe, as a function of smoothing scale, summarizing the results of the first three panels. The filled 
circles with roman numerals on the color coded line represent the radial average of the corresponding lines in panels 
A, B, and C. It is interesting to note that NGC 4254 is not unusual compared to the Local Group galaxies in terms 
of its multi-wavelength diffuse components. As expected, the FUV map has the highest diffuse fraction compared to 
the MIR 24 (im and Ha maps. The /fuv has slowest variation (~ 40 — 75%) with the smoothing scale. The /dig, 
on the other hand, falls sharply with larger smoothing scales. The amount of diffuse fraction in the 24 jam map is 
intermediate between those in the Ha and the FUV images. 

At every filter scale, we combine the DE subtracted FUV, Ha, and 24 /im emission maps to construct the desired 
Ssfr maps. Finally, a small global value is subtracted from the images. As long as most of the pixels sample the 
background, sky background subtraction can be considered an extreme case of unsharp masking where the median 
filter size is the same as the entire image. In this case, DE estimated from all three maps is negligibly small and this 
what we call as the zero diffuse fraction (/de = 0) in the remaining sections. For each Esfr map as discussed in 
Appendix [B] we use the £h 2 map derived from the combination of single-dish and interferometer data. 

To remove the diffuse extended component from the CO distribution we employ the same techniques we have applied 
to the SFR tracers. To simplify the range of parameter space to explore we use filtering kernels of the same size in 
both the Ssfr and the Sh 2 maps. Although this is not necessarily correct, we take it as an illustration of the effects 
of removing a diffuse contribution in both axes. 

F.l. Image PSF and Unresolved H II Regions 

To determine the fraction of DE in the disk of NGC 4254 it is extremely important to verify that the DE that we 
detect is a distinct source of emission in various SFR tracer maps, and that, it is not simply the emission spreads out 
from the star-forming regions. If, for example, the extended tail of the PSF contains a large fraction of the source flux, 
light from star-forming H II regions might scatter over a considerable area. This will bias the estimate of DE. 
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Figure 14. Azimuthally averaged radial profiles of diffuse fractions. Panels A, B, and C show respectively, the /fuVj /dig> an d /mir- 
In each panel, lines with Roman numerals correspond to the DE maps constructed from six different filters as given in Table [2] Panel D 
shows the average diffuse fraction (Ide) as a function of filter scale. The filled circles with roman numerals in panel D represent the radial 
averages of the corresponding lines in panels A, B, and C. The figure shows clearly that /oe increases radially outward with a trend that is 
independent of wavelength or filter scale. The profiles flatten out at the edge of the disk (beyond ~10 kpc). For smaller smoothing scales 
both /dig an d /mir vary as much as a factor of three along the disk. For longer smoothing scales the variation is about a factor of two for 
all /de- The FUV map has the highest diffuse fraction compared to the other two maps (see Table 2). Measurements for only six kernel 
sizes are shown for the ease of demonstration. The horizontal axes in panels A, B, and C span the optical radius, R25 ~12.5 kpc. 

To understand the nature of the bias we need to have clear picture about the shapes of the PSFs in the SFR tracer 
images. The PSF of these images show distinctive characteristics. For example, the PSF of Imir image has a complex 
pattern showing first and second Airy rings with radially extending artifacts. The linear scale of the second Airy ring 
from the center of the PSF is ~ 20". While approximately 85% of the total source flux is contained within the central 
peak of the PSF which has a Gaussian shape with FWHM of 6", more than 99% flux is contained with an aperture of 
diameter 40" (see Table 1 of Engelbracht et al. 2007). The PSF of GALEX NUV channel varies along the field-of-view 
from a symmetric 2D Gaussian profile to extended structure further from the center. The FWHM of the PSF contains 
more than 80% of the source flux and 95% of the total flux is contained within an aperture size of 40" in diameter 
(see Fig. 12 of Morrisset et al. 2007). The FWHM of Ihq map is fairly consistent with a 2D symmetric Gaussian. For 
this map, we find two isolated field stars far from the source and construct their light profiles. We find that 100% of 
the stellar light is confined within a region of radius 15". 

The smallest filter scale used in our study is 75". By construction, the median filter will remove any sub-structure 
of the map whose linear dimension is equal to or smaller than half of the filter scale, i.e. 38" in this case. For larger 
filter scales sub-structures of even more larger dimension will be removed. This implies that for MIR 24 /im map, at 
75" filter scale, the estimate of DE would have ~ 15% of the total emission that would come from the star-forming 
regions and not from the truly diffuse component un-associated with the region. For GALEX map one would expect 
a similar amount of contribution from the star-forming regions. For optical map it is insignificant. Note that with 
progressively larger filter scale the contribution of scattered light from H II regions becomes negligible. 

Studies of M 31 and Magellanic Clouds and Local group show that much of the diffuse materials are resolved into 
bubbles, filaments, loops, and shell like structures extending from few hundreds pc up to 1 kpc in the disk. These 
structures not only surround OB associations but also show no associations with stellar components. (Walterbos & 
Braun 1994; Kennicutt et al. 1995). In galaxies beyond Local groups the extent of most of these features along with 
faint un-resolved H II regions will fall below the resolution and hence will smear out. Therefore, any emission outside 
the bright, resolved H II regions is classified as diffuse (Kennicutt et al. 1995). Comparing luminosity functions of H II 
regions and DE Hoopes et al. (1996) showed that the contribution from the un-resolved H II regions is around 5%. 
The major contributing factor in the uncertainty of diffuse fraction estimation is the image PSF which, for this study, 
varies with the filter scale. For this analysis, however, we assign 15% uncertainty in the estimate of diffuse fraction at 
all filter scales. 



